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A composite bolometer detector system operating at 
0.37“K has been developed for making 1 mm continuum observa- 
tions using the 5 m Hale telescope on Mount Palomar. The 
cryostat employs adsorption pumping of liquid ^He to attain 
this low temperature. The electrical NEP of the bolometer 
is 1 X 10“^^ W/.'liz’. On the telescope in the best weather 
the system can detect a flux level of 5.3 Jy with a signal 
to noise ratio of 1:1 in one second of integration. 

The ^He detector system has been used to make 1 mm con- 
tinuum observations of extragalactic thermal sources. Eight 
spiral galaxies, NGC 1063, and M 82 were observed. Of this 
sample, M 82, NGC 253, NGC 4631, and NGC 6946 were detected 
with a statistical significance of better than 2a. The 3 g 
upper limits of the other galaixies were approximately 1 Jy. 
For the galaxies in the sample that have been datected at far 
infrared wavelengths (M 82, NGC 253, NGC 6946, NGC 1063, and 
IC 342) , the 1 mm points are found to lie on a modified 
black body curve 3' = 3(v,T)v^. The 1 mm results i.ndicate 
that p = 1.5 for NGC 253, M 82, and NGC 6946, while for NGC 
1068 ? ^2.0. The color temperatures derived from the far 
infrared cbseir/ations lie in the 30-40°K range. Similar 


parameters result from observations of molecular clouds in 
our galaxy. Therefore/ the observations are consistent with 
the 1 non emission from the galaxies originating in a molecular 
cloud environment. This would imply that there is not a 
'^ceat deal of 10>20"K interstellar dust in these galaxies. 

There seems to be no correlation between the 1 mm flux 
and the total HI emission from the galaxies. There is per- 
haps a correlation between the level of synchrotron emission 
and the 1 mm flux from the galaxies. 

The 1 ram flux from three points along the galactic 
plane of M 82 was measured. Only emission frcxn the nucleus 
was detected/ indicating that the dust distribution is 
strongly peaked at the location of the nucleus. This sharp 
peak is also seen in CO. 

Estimates of the gas to dust ratio for the observed 

galaxies range from less than 31:1 to greater thcui 200:1. 

These estimates assume that [CO]/[H 2 ] = 4 x 10”^ and Q/ap = 

—2 2 

5 X 10 cm /g for the dust. 


ONE MILLIMETER CONTINUUM OBSERVATIONS OF EXTRAGALACTIC 


THERMAL SOURCES 


A Thesis 

Presented to the Faculty of the Graduate School 
of Cornell University 
in Partial Fulfillment for the Degree of 
Doctor of Philosophy 


by 

Thomas Leonard Roellig 


January, 


1981 


Biographical Sketch 


Thomas Leoneurd Roellig was born on June 5, 1953, in 
Ann Arbor, Michigan. He attended public schools in Michi- 
gan cmd London, England, and received a B.S. in Physics 
from Stanford University in 1975. On August 7, 1976, he 
married Marianne Petersen. During his stay at Cornell 
University he was a teaching assistant for two years and a 
research assistant for three. In 1979 he received a Master's 
degree in Applied Physics from Cornell University. 


XI 


ACKNOWLEDGMENTS 


! 

■ i 

I 

I 

, * 

I 

I 

First and most importantly, I would like to acknowledge 
the assistance, instruction, and support of my research advis- 
or, Jim Houck. Jim's wide range of knowledge was an invaluable 
asset in my research, from the design and construction of the 
equipment to the interpretation of the results. Jim was also 
very patient and understanding, qualities that I especially 
appreciated during the observing field trips. 

I would also especially like to acknowledge the assis- 
tance of my opposite number at Cal Tech, Dave Ennis. Dave 
assisted in madcing most of the observations reported here 
and spent a great deal of his own time helping me on my vis- 
its to Cal Tech. 

The night assistants at the Mount Palomar and Siding 
Springs observatories were always willing to expend a little 
extra personal effort to make sure the observations were a 
success. I would also like to thank Gerry Neugebauer for 
the temporary use of his laboratories at Cal Tech. 

During the field trip to Australia, Dave Jauncey was 
very helpful in getting everything organised and running 
smoothly . 

George Gull and Gerry Stasavage helped in the construc- 
tion of the instrumentation. 

I would like to thank Martin Harwit and Ed Salpeter for 
serving on my committee as well as for their useful sugges- 
tions. Bill Forrest and Gary Melnick should also be 

iii 

i 


acknowledged for their discussions. 

Barbara Boettcher and Mary Wright were responsible for 
the diagrams amd typing# respectively# of this thesis. 

Last# but certainly not least# I would like to thank 
my friends amd family# and especially my wife# Marianne# 
for their support auid encouragement. 

This work was supported by National Aeronautics and 
Space Administration grant NSG-7324. 


Table of Contents 


Chapter Page 

1 Introduction 1 

2 Observing 1 mm Radiation 5 

Detectors 5 

The atmosphere 5 

Telescopes and observing sites 9 

3 Instrumentation 12 

Cryostat 12 

Bolometer 20 

Filters 21 

Photometer 21 

Telescopes 22 

Electronics and data system 23 

4 Observing Procedures and Data Reduction . . 30 

Observing procedure 30 

Millimeter flux standards 37 

Data reduction 39 

5 Observations of Extragalactic Thermal 

Sources 41 

Introduction ..... 41 

Observations 42 

Millimeter and submiliimeter spectra . . 43 

Dust masses 65 

Flux correlations 70 


V 


Table of Contents (cont'd) 


Chapter Page 

Individual objects 79 

M 82 79 

NGC 253 85 

NGC 6946 87 

NGC 4631 88 

IC 342 90 

NGC 1068 91 

upper limits 99 

The Milky Way galaxy 100 

Summary 102 

Appendix A The Bolometer 104 

Theory and performance 104 

Total system noise Ill 

Bolometer construction 119 

Appendix B Filters and Field Optics 129 

Appendix C Atmospheric Corrections 142 

References 151 


List of Tables 


Table Page 

3-1 Cryostat Performance 20 

3-2 Bolometer Performance 21 

3- 3 Preamplifier Performance 26 

4- 1 Calibration Standards 38 

5- 1 Nuclear Positions 44 

5-2 Observations 45 

5-3 Isothermal Oust Model Parcuneters 59 

5-4 Two Temperature Dust Model 65 

5-5 Values of Q/ap at 1 mm for Some 

Suggested Grain Materials 69 

5-6 Oust Mass Column Densities 71 

5-7 Galactic Distances 72 

5-8 Total Dust Masses 73 

5-9 Upper Limits 100 

A-1 Symbol Definitions 107 

A-2 Dewar Noise at 10 Hz, Measured 

at the Bolometer 113 

A-3 Measured Bolometer Electrical Performance. . 113 

A-4 Heat Capacity of Bolometer Elements 122 

B-1 Bolometer Optical Performance 141 

vii 


List of Illustrations 


Figure Page 

2- 1 Atmospheric transmission and dewar 

response * . . . . 7 

3- 1 ^He dewar 14 

3-2 Pumping system. 16 

3-3 Beam profile 24 

3-4 5 m data system 27 

3- 5 4 m data system 29 

4- 1 Beam nomenclature 31 

5- 1 Galactic spectra 47 

5-2 Flux derivation nomenclature 67 

5-3 HI -1mm plots 76 

5-4 Synchroton - 1 mm plot 78 

5-5 Scans along plane of M 82 82 

A-1 Schematic diagram of a bolometer 105 

A-2 ^He bolometer signal and noise 114 

A-3 Noise spectrum 116 

A-4 ^He bolometer , 121 

A-5 Bolometer soldering 126 

B-1 Filters and field optics used in dewar. . . . 130 

B-2 Diffraction transmission curves . 134 

B-3 Dewar response -grating monochromater .... 139 

B-4 Dewar response - Fourier transform 

spectrometer 140 


viii 


List of Illustrations (cont'd) 


Figure Page 

C“1 Q curves 144 

C~2 Corrected flux of 7027 145 

C-3 Effective wavelength of system. 149 


Chapter 1 


Introduction 

The field of millimeter continuum astronomy is still 
in its infancy. Observations are so difficult that the 
technology needed to make useful astronomical measurements 
has only recently been developed. The difficulty of making 
observations is due to three factors. 

First of all, the atmosphere is essentially opaque to 
radiation from 20 ym to 750 um as far as ground based 
observations are concerned. Beyond 750 ym the atmosphere 
becomes more transparent, but the atmosphere still attenu- 
ates and introduces noise into the celestial radiation. 

The second difficulty in making 1 mm observations is 
caused by the state of the art of detector technology. 
Detectors in the near and mid infrared regions of the spec- 
trum (1 ym-30 ym) have a great deal of use in military and 
industrial applications. As a result, development of these 
detectors has been well funded and the detector technology 
is relatively advanced. Typically, -nfrared astronomers 
working in this wavelength region are able to purchase 
excellent detectors from industrial sources. The detectors 
used in 1 mm astronomy are generally fabricated by the 
astronomers themselves and have not been subject to the 
intensive development of the photoconductor and photovol- 
taic devices used in shorter wavelength work. As a result 
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of this limitation, 1 mm astronomy is as much a field of 
detector technology as it is observational astronomy. 
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Lastly, most astronomical sources are very faint at a 
wavelength of 1 mm. This wavelength region is long enough 
to be well into the Rayleigh-Jeans tail of thermal sources 
(T > 5®K) and yet too short for strong synchrotron emission. 
In spite of these difficulties, 1 mm observations are im- 
portant since they Ccin provide information about astronomi- 
cal objects that is unavailable at other wavelengths. 

Astronomical objects that are observed at 1 ram can be 
put into three main categories. 

a) Thermal sources. 1 mm thermal sources include 
solar system bodies, spiral galaxies, molecular clouds. 

Class II Seyfert galaxies, and planetary nebulae. With the 
exception of the planetary nebulae all of these sources 
emit 1 mm radiation either from the surfaces of the objects 
^hemselves or from heated dust grains. The 1 mm emission 
from planetary nebulae is due to thermal bremstrahlung . 

The brightest sources in the 1 mm sky, the sun, Jupiter, 
and Saturn, are all thermal sources. 1 mm observations of 
the solar system objects can provide information about 
their atmospheres and the composition of their surfaces. 
Measurements of the 1 mm flux from planeteiry nebulae can 
check for the presence of ultra compact components. As 
will be seen in more detail later, 1 mm measurements and 
maps of dusty objects can provide very useful information 
about the temperature, composition, and distribution of the 
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dust. 

b) Non-thermal sources. Observed non-thermal sources 
have been either radio galaxies, QSO's, Class I Seyferts, 
or supernovae renmants. In general, the 1 min flux from 
these objects can be predicted by extrapolating from the 
radio spectrum (Ennis, 1980) . Measurements at 1 mm are use- 
ful because they can set limits on the energy density and 
magnetic fields in these objects. 

c) Cosmic baclcground. The 3*K cosmic background has 
roughly a blackbody spectrum. Thus, its peak is near a 
wavelength of 1 mm. Unfortunately, the severe atmospheric 
noise generally negates the value of photometric observa- 
tions at this wavelength, since accurate calibration is 
impossible. As a result, most cosmic background work has 
been in the context of investigations of the spectrum. 1 mm 
observations are also useful for special projects such as 
searching for small scale inhomogeneities in the background 
radiation or observing the Sunyavev-Zeldovich effect in 
clusters of galaxies (Gould and Rephaeli, 1978). 

This work reports the results of 1 mm observations of 
extragalactic thermal sources. The methods of making 1 mm 
observations are described in Chapter 2. The instrumenta- 
tion used to make the observations is described in Chapter 
3, with more details on its theory of operation and con- 
struction given in Appendices A and B. The observing proce- 
dure followed is given in Chapter 4. More details on the 
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data reduction, especially the atmospheric correction pro- 
cedures, are given in Appendix C. The results of the 
observations, together with their interpretation, is given 
in Chapter 5. 

This work was done in collaboration with members of the 
infrared group at Cal Tech. In this collaboration the 
Cornell group provided the detector and filters, while the 
Cal Tech workers provided the photometer. We jointly 
obtained the telescope time on the 5 m Hale telescope. The 
observations themselves were made by joint research teams 
composed of members of both groups. The observations and 
results reported here constitute part of the author's con- 
tributions to the program. The other results of this col- 
laboration appear iix the open literature. 


Chapter 2 


Observing 1 mm Radiation 


Detectors 

Two different types of detectors are presently in use 
in making 1 mm astronomical observations. The first type 
of detector is a superheterodyne system similar in princi- 
ple to those used in radio astronomy. In millimeter work 
an InSb crystal usually functions as the mixer. The sensi- 
tivity of these devices is limited by the power of the 
local oscillator. Since high harmonics of a Klystron 
generator are used to supply the local oscillator power, 
the sensitivity of these devices is relatively low, with 
system temperatures of a few thousand degrees. Since the 
bandwidth of these devices is quite small, they have been 
used chiefly in molecular line work. 

The second type of 1 mm detector is the bolometer. 

These devices are incoherent detectors with very large 
bandwidths. As a result, they are more suitable for con- 
tinuvim work thcin cire the superheterodyne receivers . The 
detector used in this study was a bolometer. The theory 
of operation for these detectors and the methods of their 
construction are treated in Appendix A. 

The Atmosphere 

One of the biggest problems in making 1 mm observations 
is caused by the atmosphere. As can be seen in Figure 2-1, 
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the atmosphere is essentially opaque at ground based observ- 
ing sites from 20 ym to 750 ym. Most of the atmospheric 
absorption is due to water vapor, with smaller amounts of 
absorption due to oxygen, ozone, and carbon dioxide (Traub 
and Stier, 1976) . 

In reducing the 1 mm results, the first step is to 
remove obvious noise spikes from the data. The data is then 
corrected for the effects of the atmosphere. In practice, 
this means usin(a an atmospheric model to find the effective 
attenuati oir. through the detector wavelength passband. 

Since the dominant source of atmospheric opacity is due to 
water at these wavelengths, good atmospheric correction 
procedures need accurate water vapor measurements. 

One way of measuring the water vapor content of the 
atmosphere is to measure the dew point during the observa- 
tions. This technique has been employed by Rowan-Robinson 
et al. (1975) in reducing their 1 mm data. Because this 
method only measures ground level water vapor, non-uniform 
atmospheric conditions cause inaccuracies in this procedure. 

A more accurate method is described in Elias et al. 
(1978) and was used for the 1 mm results presented here. 
Briefly, this method employs measuring the infrared flux 
from the sun both in and just outside of the 1.9 ym water 
band. This measurement is made using a calibrated Westphal 
meter (Westphal, 1974) . This technique has the advantage 
that it measures the integrated water vapor content through 
the whole atmosphere. Unfortunately, with the present 
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Figure 2-1; 


a) Atmospheric transmission from 1 urn to 
1000 urn, with 1.2 mm of precipitable water 
vapor (Traub and Stier, 1975). 

b) Smoothed passband of detector and dewar 
used in the 1 mm observations reported 


here. 


100 ^ 
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equipment the method cannot be used at night. The Westphal 
meter readings can be converted into precipitable milli- 
meters of water, which can then be used as a parameter in 
an accurate atmospheric correction model. The atmospheric 
model used to reduce the data presented in this work is 
described in Appendix C. 

Telescopes and Observing Sites 

1 mm observations have been made from balloons cu>.d 
ground based observatories. Aircraft observatories have 
been used in some submillimeter observations, but have not 
been used at 1 mm. The use of high altitude balloons has 
the advantage of removing a great deal of the atmospheric 
attenuation of the celestial signals, permitting observa- 
tions at submillimeter wavelengths as well as at 1 mm. 
Balloon telescopes have the disadvantage of small size 
(generally their mirror diameters are 1 m or less) , which 
limits both the spatial resolution and the collecting area. 
They also have the disadvantage that their flights are of 
limited duration, thus long observing sessions are not pos- 
sible. As a result of these drawbacks, balloon telescopes 
have only been used at millimeter wavelengths for 3°K cos- 
mic background work, where the reduction of atmospheric 
effects is the overriding concern (Woody and Richards, 1979). 

For most of the past work at 1 ram, ground based tele- 
scopes have been used. These telescopes are either the 


10 


millimeter dishes used in molecular line radio astronomy or 
large optical telescopes. Large ground based telescopes 
can be much larger than balloon borne telescopes. Their 
diffraction limited beeun sizes are of the order of 1' , 
which allows reasonably fine scale mapping to be undertaken. 
The large apertures also make the detection of faint point 
sources possible. There are two disadvantages in using 
these telescopes. First of aJl, their use is subject to the 
weather. Since observations are generally made during the 
winter so that the water vapor levels will be low, the 
weather is frequently cloudy. The second disadvantage lies 
in the popularity of these telescopes; it is difficult to 
get observing time. 1 mm astronomers are fortunate in that 
they cam work in the daytime, so that they can share the 
telescope with a nighttime observer. 

For ground based work at 1 mm, the observing site is 
just as important as the telescope that is used. The ber^t 
sites are at locations where the water vapor content of 
the atmosphere is very low. In practice, this means that 
the observatories are located at high altitudes. For 
example, as far as sensitivity is concerned, the 5 m Hale 
telescope on Mount Palomar (alt. 1706 m) is much better than 
the 10 m dish at Owen's Valley (alt. 1216 m) . This is due 
to a smaller amount of water vapor over the Hale telescope, 
which more than balances out the factor of four in the 
collecting area. The ideal site for 1 mm ground based 
astronomy is atop Mauna Kea in Hawaii, since this location 
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has very little water vapor. Plans exist to build a 10 m 
millimeter (Leighton) dish at this location in the near 
future. Longer range plans include building a 25 m tele<- 
scope for Mauna Kea. The observations presented in this 
work were made at either the 5 m Hale telescope or at the 
4 m Anglo-Australian telescope (alt. 1165 m) . 


Chapter 3 


Instrumentation 

The 1 mm observations reported here were made using a 

3 

composite bolometer operating at pumped He temperatures. 

The reasons for going to such low temperatures are detailed 
in Appendix A. This chapter is concerned mainly with describ 
ing the cryostat/ photometer, telescopes, and electronics of 
the instrumentation. 

Cryostat 

•3 

For astronomical applications, using pumped “'He cryo- 
stats is the most convenient way of reaching temperatures of 
less than l.C’K. Such cryostats are compact, relatively 
simple, and stable, facilitating their use on telescopes. 

Pumped ^He reaches lower temperatures than pumped ^He for 

^ 4 

three reasons; a) -^He has a higher vapor pressure than He 

at all temperatures, thus reaching lower temperatures for a 

given pumping pressure, b) ^He is a fermion, so there is no 

superfluid film to transport heat for T > 10 m“K, and c) the 

lack of a superfluid film means that the pumping lines can 

be constructed of large dieuneter tubing without the small 

orifices needed in pumped ^He dewars. Since ^He is an 

expensive gas (- $160/liter STP) , the pumping mechanism 

employed must be closed so that the ^He can be recovered. 

For our purposes, an additional constraint in the dewar 

design was set by size, since the dewar was intended to be 
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compatible with an existing photometer on the 5 m Hale 
telescope. 

The basic dewar configuration is shown in Figure 3-1 
and consists of a liquid N2 bath and radiations shield, a 
liquid He bath and work surface, and a vacuum jacket. The 
maximum height of the dewar is only 16". 

The ^He pumping system is shown in Figure 3-2. An 
adsorption pump with a mechanical heat switch was used to 
provide a fast, closed, and vibrationless pump, with the 
exception of the condenser, the pumping lines running from 
the liquid chamber to zhe adsorption pump are all of thin 
wall 3/16" stainless steel tubing, giving a pumping speed at 
4*K greater than 1 1/2" tubing at room temperature. The con- 
denser is constructed of OFHC copper and is thermally 
clamped to the ^He work surface. In order to give a good 
thermal contact between the work surface and the condenser, 
the surfaces that were clamped together were gold plated. 

The gas flows through the condenser through seven 1.6 ram 
holes. This impedance, although slowing the pumping speed 
slightly, increases the cold surface available for condensa- 
tion and stops oscillations in the ^He pumping line. The 
adsorption pump consists of a 2" dieuneter OFHC copper ball 
filled with activated charcoal. The ball has numerous inter- 
nal fins made of OFHC copper so that ai^y piece of charcoal 
is in close proximity to a cold surface. A 1/8" thin wall 
stainless steel tube runs from the ball to a refrigeration 
gauge and various plumbing fittings on the top of the dewar. 



Figure 3-1: Drawing of dewar used in observations* 
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Schematic diagram of the 
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He adsorption 


pumping system. 
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The refrigeration gauge is used for monitoring the ^He pres- 
sure, while the fittings are used to charge the dewar with 
■^He. The heat switch is mechanical, with a geared drive 
that strongly presses two gold plated OFHC copper contacts 
together. One of these contacts is soldered to the ^He work 
surface while the other is soldered to a thick braided 
copper grounding strap connected to the ball containing the 
activated charcoal. The heat switch shaft leading into the 
dewar from the outside world has an 01dh^un coupling. This 
allows the shaft to be disconnected to reduce the thermal 
heat load into the dewar. 

The operation of the cryostat is as follows. 

1) The ^He bath is topped off and the heat switch is 
opened, breaking the thermal path between the ball and the 

4 

He work surface. 

2) A special fitting is attached to the ^He fill tube 
and the ‘^He bath is pumped to as low a temperature as possi- 
ble. At the same time, current is run through a 100 ohm 
heating resistor on the ball. A 15 volt input puts 1 watt 
of heat into the charcoal. (The missing power is dissipated 
in the wires inside the dewar, which have a resistance of 
50i^.) The heater is turned off about 25 minutes later when 
the charcoal reaches approximately 50 °K. This temperature 
is high enough to insure that all of the ^He is driven from 

4 

the charcoal. As soon as the He bath drops below the con- 
densation temperature, the ^He will start condensing in the 
condenser. At 1 atm the condensation temperature of ^He is 
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3.2®K. The critical temperature is only 3.35®K, so that the 
^He bath must be pumped to a lower temperature than this 
before condensation occurs. Gravity causes the ^He droplets 
to then run down the pumping line and into the liquid chamber. 
Eventually, the ^He bath is as cold as it will get (about 
1.5“K in the present dewar configuration) . At this time 
approximately 90% of the ^He is in liquid form, filling a 
volume of 1.4 cm^ in the liquid chamber. 

3) The heat switch is then closed, which starts cool- 
ing the adsorbant. At the same time the pump is disconnected 
from the ^He bath, allowing the ^He bath temperature to rise 
to 4.2“K. As the ball temperature drops below 15*K the ^He 
will start to readsorb onto’ the charcoal. Within 20 minutes 
the temperature in the liquid ^He chamber will have dropped 
to the operating temperature. This temperature has been 
determined to be 0.37®K by a calibrated Ge thermometer. The 
loss in cooling the ^He to this temperature is only about 15% 
of the original volume of liquid. The ^He itself has the 
greatest heat capacity and its self-cooling is responsible 
for most of this loss. 

4) Topping up the ^He bath completes the start-up 
procedure. 

The heat load on the ^He is low enough so that the hold 
time at 0.37'’K is approximately 15 hours. The complete list 
of the cryostat performance specifications is given in Table 


3-1. 
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Table 3-1 

Cryostat Performance 


Liquid N 2 

hold time 

12 hours 

Liquid He 

hold time 

14 hours 

Pumped ''He hold time 

15 hours 

Operating 

temperature 

0.37»K 

Heat load 

into the ^He 

13 uW 

Thermal conductivity of heat switch 

45 mW/®K 

Volume of 

3 

He pumping system 

120 cm^ 

Amount of 

^He used in dewar 

1 liter STP 

Volume of 

^He bath 

0.87 liters 

Volume of 

liquid N 2 bath 

0.80 liters 


Bolometer 

A full discussion of the bolometer, including the 
theory of operation, performance, and construction, is given 
in Appendix A. Briefly, the bolometer is a composite design 
with a semiconductor thermistor element. A listing of its 
performance is given in Table 3-2. On the 5 m Hale tele- 
scope in the best weather the bolometer can detect 5.8 Jy 
with a signal :noise ratio of one in 1 second of integration 


time 
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Table 3-2 

Bolometer Performance 


Electrical NEP 
Optical NEP at 10 Hz, f/35 
Thermal time constant 
Operating temperature 


1.5 X 10“^^ W//Hz 
beam 5 x 10”^^ W//Hz 
14 ys 
0.37“K 


Filters 

The filters used for the observations are covered in de- 
tail in Appendix B. We will only note here that the half 
power points of the dewar transmission passband are at 750 ym 
and 1200 ym. 

Photometer 

The bolometer dewar is mounted on a chopping photometer 
at the prime focus of the telescope. In its original con- 
figuration this photometer has been described in Elias et al. 
(1978). Briefly, it is a reimaging system with a wobbling 
tertiary mirror that performs the chopping. Recently, the 
photometer has been modified at Cornell for a larger chopper 
throw and for remote control operation. 

The photometer reimages the beam from the telescope into 
the dewar. At the same time, the beam shape also changes 
from the telescope's f/3.3 to f/4.0. The photometer servo 
drives the tertiary mirror to follow an input waveform 
provided by a standard function generator. The photometer 
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can chop as far as 4' and at a rate of 20 Hz before the mir- 
ror encounters difficulty in following a square wave input. 
The photometer orientation can be changed to chop in any 
direction. At the 5 m Hale telescope the drive electronics 
only permit the telescope to nod automatically in the N-S 
or E-W directions. 

A motor driven 45® flat mirror is located inside the 
photometer top. With the mirror in the beam the light is 
reflected into a television camera. In the viewing mode 
accurate telescope pointing can be obtained from observations 
of bright stars. The television camera has enough sensitiv- 
ity and contrast so that stars down to the 4th magnitude can 
be observed in daylight against a deep blue sky using the 
5 m Hale telescope. This is about as well as the human eye 
can do. At night the television system limiting magnitude 
is roughly 8.5. The human eye does much better than this! 
Therefore, for certain applications where pointing must be 
determined from dim stars at night, the television Ccimera 
can be removed and an observer can ride in the prime focus 
cage and do the pointing. In practice, the television sys- 
tem has been sufficient for all of our observations. When 
the mirror is pulled out of the beam, the radiation enters 
the bolometer dewar. The mirror position may be changed from 
the observatory data room. 

Telescopes 


The observations were made at either the 5 m Hale 
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telescope on Mount Palomar or the 4 m Anglo-Australian tele- 
scope at Siding Springs. The same photometer was used on 
both telescopes/ since the beam shapes at prime focus were al- 
most identical. The dewar incorporated a field stop of 1/4" 
diameter. With the f/4 be6un on the 5 m telescope the geomet- 
rical beeun size was 66" on the sky. Due to diffraction ef- 
fects, the beeun shape was not perfectly square, but could be 
approximated by a Gaussian core with enhanced wings (Figure 
3-3). The FWHM of the measured beam was 55". The Anglo- 
Australian telescope was similar to the Hale telescope, except 
that the primary mirror diameter is 3.8 m. As a result, the 
geometrical beam size on the Anglo-Australian telescope was 
88" on the sky. Lack of a suitably bright source precluded 
measuring the beam shape on the 4 m telescope; it was assumed 
to have been the same as on the 5 m telescope. This would 
correspond to a F'»iJHM beam size of 73". The observing proce- 
dures used on these telescopes are given in Chapter 4. 

Electronics and Data System 

The voltage signals from the bolometer are first am- 
plified by a preamp mounted on the dewar. The performance 
of this amplifier is shown in Table 3-3. The amplifier 
was designed to have the lowest noise possible. The JFET 
used in the front end was specially selected for its low- 
noise characteristics. To further reduce the noise, the 
preamp is powered by two 13.6 volt mercury batteries. The 
capacity of the batteries is large enough so that the pre- 
amp can be operated for over 40 hours before battery changes 



Figture 3-3 : Beam shape of the 1 mm system on the 5 m 

Hale telescope, obtained by scanning the 
planet Saturn* 
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Table 3-3 

Preamplifier Performance 


Voltage gain at 10 Hz 
Voltage noise at 10 Hz 
Current noise at 10 Hz 


858 

4.1 nV//Hz 
3.4 fA/^Hz 


become necessary. The bolometer bias voltage source is also 
contained in the preamp box and is powered by these seune 
batteries. Proper shielding is very important in eliminating 
pickup from nearby electronics. In order to provide this 
shielding the preamp was bolted directly to the vacuum jacket 
of the bolometer dewar. In addition, the wires connecting 
the front end of the amplifier to the bolometer and load 
resistor were made as short as possible and were doubly 
shielded. 

The output of the preamp entered either of two data sys- 
tems (shown schematically in Figures 3-4 and 3-5) depending 
on the telescope in use. At the Hale telescope, integration 
was performed by counting the cycles out of a voltage to fre- 
quency converter. At the Anglo-Australian telescope the 
in-house computer system allowed us to perform the integra- 
tions in the data acquisition computer after sampling the 
output of the lock-in amplifier. 


Figxaxe 3- 
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: Schematic 


diagram of the data system used 


on the 5 m Hale telescope. 























Chapter 4 


Observing Procedure and Data Reduction 
Observing Procedure 

During the 1 mm observations special observing proce- 
dures were followed to insure low noise accurate fl\ix meas- 
urements. Because the atmosphere causes a great deal of 
trouble in ..iaki-'g 1 mm astronomical observations, most of the 
observing procedures are therefore employed to minimize the 
effects of the atmosphere. 

The observed millimeter flux from any object aside from 
the sun is much less than the emission from the atmosphere 
and telescope. As a result, differential measurement tech- 
niques that allow partial cancellation of the atmospheric 
emission are used. When the wobbling tertiary mirror in 
the photometer is driven by a square wave signal, it has the 
effect of rapidly alternating the telescope position between 
two spots in the sky. This is shown in Figure 4-1. When 
the telescope is positioned so that the object being observed 
is in the left beam, the telescope is defined to be in the 
signal position. Conversely, when the telescope is posi- 
tioned 30 that the object is in the right beam, the tele- 
scope is defined to be in the reference position. If Fj^g 
and Fl 3 are the fluxes entering the bolometer dewar from the 
left and right beams, respectively, we find that the AC sig- 
nal when the telescope is in the signal position is 
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Figure 4-1: Diagram of the beam position-.^ and nomenclature 

used in differential flux measurements. In the 
signal position the object is in the left beam 
as the wobbling tertiary mirror chops in a 
square wave pattern at 11 Hz. To move to the 
reference telescope position the telescope 
is nodded the distance of the chopper throw. 

The telescope nods back and forth every 20 


seconds. 
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Ss = <4.1) 

where R is the dewar responsivity in volts per watt. Breaking 
down the total fluxes into their various components we have 

Sg = R((Fo + Pt)i,B - (Pa + f't'RB) (<-2) 

where Fq is the flux from the object, is the flux from 
atmospheric emission, and F^ is the flux from emission from 
the warm telescope and photometer. Similarly., the AC signal 
when the telescope is in the reference position is 

Sr = " FrbI (4.3) 

or 

Sr = RKPa + Pt)LB - (fa + + Vre’ (‘'•‘" 

The electronics of the data system synchronously demod- 
ulate and integrate the AC signals. As long as the atmos- 
pheric emission is the same in both beams (a good approxi- 
mation except when there are clouds or an approaching frontal 
system) the signal from the object will then be given by 

So = 1^0 = f I /Sg - / Sr] (4*5) 

This procedure therefore cancels out almost all of the 
effects of the telescope and atmospheric emission. 

The telescope is nodded between the signal and reference 
positions every 20 seconds. This allows two 10 second inte- 
grations in each position. When the weather is only marginal, 
it is sometimes advantageous to nod the telescope more 
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frequently, perhaps every three seconds. This tends to re- 
move some of the excess sky noise. The disadvantage with 
shortening the integration time is that the time it takes 
the telescope to move and settle aiiter every nod is time 
lost for integration purposes. A standard run for a very 
bright object such as a planet is two pairs, or SRRS, where 
S and R refer to integrations with the telescope in the sig- 
nal and reference positions, respectively. For longer 
integration runs a SRRSSR...RS pattern was followed. This 
had the effect of cancelling out the effects of a monotonic 
linear change in the weather conditions. The standard ob- 
serving block for fainter objects was composed of 16 pairs 
of 10 second integrations. 

The details of the observing procedure are giv^n below. 

1) Initial set-up ; The first task undertaken in any 
observing trip is to allign the millimeter and optical beeims. 
To do this it is necessary to observe an object that is 
bright at both 1 mm and optical wavelengths. The ideal 
objects are the planets Venus, Mars, Jupiter, and Saturn, 
since they are not only bright in both wavelength regions 
but are also almost point sources for our 55" becun. This 
last feature means that the beam center can be easily deter- 
mined. 

With the 45® flat in the photometer head out of the 
beam, the 1 mm beam center is found by observing the output 
of the lock-in aunplifier while moving the telescope. The 
strip chart connected to the lock-in amplifier output makes 


35 


this task easier. When the 1 mm signal peak is found, the 
telescope is then moved to this position and the 45" mirror 
is inserted into the beam. The light from the object is now 
reflected into the television camera. An adjustable reticle 
located in the photometer head is then moved so that the 
cross hairs are centered on the left-hand image. For the 
rest of the observing run any object centered in the cross 
hairs will also be centered in the signal 1 mm beam. 

The 45“ mirror is next moved out of the beam and the 
phase control on the synchronous demodulator and the tele- 
scope focus are set to provide a maximum 1 mm signal. The 
telescope is then nodded into the reference becim. If the 
signal becomes negative, all is well. If the signal merely 
drops to zero, the object was centered in the reference 
rather than the signal beam. The phase on the synchronous 
demodulator should be changed by 180“ and the signal becim 
located by moving the telescope and inspecting the strip 
chart as before. 

The chopper frequency is also set at this time? the fre- 
quency was chosen to give the highest signal to noise. In 
practice, using the Cornell ^He bolometer dewar and with 
average atmospheric noise, this frequency was around 10 Hz. 

In order to eliminate any sort of beating with the 60 Hz 
mains the frequency was set to some odd number, such as 11 
Hz. 

2) Observing a standard ; Before any set of observa- 
tions it is necessary to observis some millimeter standard so 
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that the atmospheric transmission and instrumental sensitiv- 
ity can be determined. The primary standards for this purpose 
are the planets# since their 1 mm brightness temperatures 
have been measured by Werner et al. (1978). Other secondary 
standards ceui be used when planetary observations are impos- 
sible. These standards will be discussed in more detail 
shortly. In general, measurements of a standard source are 
made approximately every 45 minutes, or even more frequently 
when the atmosphere appears to be unsettled. In addition, 
readings of the water vapor using the Westphal meter are 
made once every hour that the sun is up. 

3) Observations of the object; Determining accurate 
telescope pointing is the first step in observing a milli- 
meter object. The easiest way to do this is to offset point 
from a bright star that is close to the position of the 
millimeter object. With the 45® mirror in the beam the left- 
hand image of the star is lined up in the reticle cross 
hairs. The telescope control computer is given commands 
that update the position readout to match the precessed 
catalogued position of the bright star being observed. After 
these corrections hc've been made, the telescope is set on 
the precessed position of the millimeter object. Finally, 
after the telescope is correctly positioned, the 45® mirror 
is moved out of the beam and integration is started. 

During any observing trip step 1 is done once (or maybe 
twice to check that it was done correctly the first time) at 
the beginning of the run. The other i-teps are done with 


each new object. If the integrations on an object will be 
particularly long, step 2 may be done a few times in the 
middle of the integration run to check for weather changes 
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Millimeter Flux Standards 

Since the atmospheric effects are so large at 1 mm, 
extraterrestrial objects are needed for absolute flux cali- 
bration standards. This poses a problem, since both the 
temperature and emissivities of astronomical objects are un- 
known. For the 1 mm results reported here we have used the 
planets Venus, Jupiter, and Saturn as our primary standards. 
The 1 mm brightness temperatures of these bodies were taken 
from Werner et al. (1978). In this paper, the reported 
temperatures were all derived from measurements of the 
planets against the planet Mars. The 1 mm brightness tem- 
perature of Mars was calculated from Mariner 9 results using 
the thermal model of Wright (1976) . Recently, there has 
been evidence that this model may be incorrect (Forrest, 
1980). Taking this into account, we have adopted the values 
shown in Table 4-1 for the 1 mm brightness temperatures of 
the planets. 

When the planets were observed, corrections had to be 
made to correct for the partial resolution of the planetary 
disk by the beam. This correction is given by the formula 


So 





(4.6) 
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Table 4-1 

Calibration Standards 


1 mm brightness 
276 
168 
145 
87 


Primary standards 

Object 

Venus 

Jupiter 

Saturn 

Uranus 

Secondary standards, 

Object 
NGC 7027 

3C84 (Dec. 1 , 1979) 
(Dec, 1 , 1979) 
(Dec. 30, 1979) 
(Feb. 29, 1980) 
(Mar. 29, 1979) 
(Mar. 29, 1979) 


derived fluxes 

1 mm Flux (Jy) 
5.8 + 0.8 

53.2 + 6.1 

41.2 + 6.7 

41.1 i 6.7 
50.9 + 5.8 

63.2 £ 7.5 
56.8 ± 8.4 


temperature (“K) 
± 31 
± 19 
± 16 
+ 11 


Primary Standard 
Saturn 
Jupiter 
Saturn 
Jupiter 
Venus 
Venus 
Saturn 


Note: The measurements of 3C84 made against Saturn were 

actually made against NGC 7027, but the calibration of NGC 
7027 came from Saturn. 
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where Sq is the observed flux from the planet, is the 
actual flux from the planet, d is the semidicimeter of the 
planet in arc seconds, and a is related to the FWHM of the 
beam shape by 


g - FWHM/2 

1/2 


(4.7) 


These standard sources are then used to calibrate other 
bright objects, such as 3C84 and NGC 7027 so that they can 
be used as secondary standards. During the period of these 
observations, 3C84 began to flare. The fluxes measured from 
3C84 during the different field trips are therefore listed 
separately in Table 4-1. The 1 mm flux measured for NGC 7027 
was in good agreement with previous values from Elias et al. 
(1978). Since the temperature of Mars varies throughout the 
year, it was also measured separately against one of the 
giant gas plcuiets on each field trip that it was used as 
a calibration source. 


Data Reduction 

As was noted above, data from faint sources was general- 
ly taken in blocks of 16 pairs of integrations. This data 
was reduced in the following way. First, the strip chart 
recorder output was carefully scanned and any obvious noise 
spikes were removed from the data base. The cleaned data 
was then reduced by dividing out the atmospheric effects 
using the Q's described in Appendix C. Typically runs of 
four blocks of 16 pairs each were sandwiched between 
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calibration observations of primary or secondary standards. 

After cleaning and correction for atmospheric effects/ 
the data from the four blocks was then averaged together/ 
and a mean and statistical error computed. The resulting 
number was then expressed as a fraction of the cleaned and 
corrected signal from the calibration runs from the standard. 
Multiplying this fraction against the 1 mm flux from the 
standard then gives the 1 mm flux from the object. 


Chapter 5 


1 mm Observations of Extragalactic Thermal Sources 
Introduction 

The extragalactic thermal sources that could be expected 
to be detected at 1 millimeter with today's technology are 
either spiral, irregular, or class II Seyfert galaxies. 
Elliptical galaxies have little gas and dust compared to 
spiral galaxies and their synchrotron emission is relatively 
weak. As a result, no elliptical galaxy is expected to be 
detectable at 1 mm using present instrumentation. 

Some spiral galaxies do have enough 1 mm flux to be 
detected with our ^He detector. For normal spiral galaxies 
the 1 mm radiation due to dust emission is stronger than the 
synchrotron or thermal bremstrahlung emission. This can be 
seen by examining the spectrum of NGC 253, a typical spiral 
galaxy, shown in Figure 5-la. 

Seyfert galaxies have been divided into two classes on 
the basis of the width of their forbidden emission lines 
(Khachikian and Weedraan, 1971, 1974). There is a funda- 
mental difference in the infrared emission of the two 
classes. The infrared radiation from the class II Seyferts 
is thought to arise predominantly from thermal emission from 
dust grains, while the infrared radiation from class I Sey- 
ferts is due to the srme non- thermal mechanism that is re- 
sponsible for the observed short wavelength radio emission 
(Neugebauer et al. , 1976). Both types of Seyfert galaxies 
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have been observed at 1 mm in the past, with only class I 
Seyferts being detected (Elias et al., 1978; Hildebrand et 
al. , 1977) . 

Observations 

This study reports observations of 8 spiral galaxies, 

M 82, euxd NGC 1068. These galaxies were selected on the basis 
of the following criteria. First of all, the galaxies were 
selected so that their nuclear regions would fill, or nearly 
fill, the 1' beam. In addition, edge on galaxies were favored 
because of the greater path length through their galactic 
disks. Finally, the galaxies that had been observed at other 
wavelengths, especially in the far infrared and in CO, were 
favored. 

The observations were carried out at the 5 m Hale tele** 
scope at Mount Palomar and the 4 ra Anglo-Australian tele- 
scope at Siding Springs. In order that the v/ater vapor con- 
tent of the atmosphere be as low as possible, the observations 
were made during local winter. The weather at the 4 m 
telescope during observations of NGC 253 was only marginal, 
with 5 mm of precipitable water vapor and a high level of 
sky noise due to an advancing frontal system. The observa- 
tions at Mount Palomer reported here were made during good 
to excellent weather conditions. 

The positions of the galactic nuclei were measured from 
Palomar Sky Survey plates, or in some cases, were taken from 
Dressel and Condon (1976) . The positions of the galactic 
nuclei are given in Table 5-1. Three positions in M 82 were 
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observed, the nuclear position and points 1' on each side of 
the nucleus along the plane of the galaxy. 

The results of the observations of the galaxies are 
shovm in Table 5-2. Two different errors are shown in this 
table. The errors associated with the fluxes are the total 
uncertainties in the fluxes, including statistical errors in 
the observations of the object and millimeter standards, 
as well as uncertainties in the 1 mm flux of the standards. 
The errors listed in column 6 are statistical errors from 
the observations of the object only, expressed as Signal/RMS 
Error. Due to uncertainties in the observations of the stan- 
dards and the uncertainty in the known fluxes of these 
standards, these errors are smaller than the total errors 
given in column 5. 

Millimeter and Submillimeter Galactic Spectra 

i) Temperatures ; By combining the 1 mm continuiim ob- 
servations with shorter wavelength data, it is possible to 
estimate the dust temperatures and gain information about 
the grain composition. This information will be needed to 
estimate the dust masses of the gal^ucies. 

As a first approximation, the submillimeter spectra 
can be fit with curves of the form 

F(v) = e{v)B(v, T) (5.1) 

where B{v, T) is the Planck blackbody function and s(v) is 


of the form 
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Table 5-1 
Nuclear Positions 


Object 
M 82 nucleus 

9 

1' NE 

9 

1' SW 

9 

NGC 253 

0 

NGC 891 

2 

NGC 925 

2 

NGC 1068 

2 

NGC 4244 

12 

NGC 4631 

12 

NGC 5907 

15 

NGC 6946 

20 

IC 342 

3 


RA (1950) Dec 
h 51 m 45.6 s 69® 
h 51 m 55.7 s 69® 
h 51 m 35.5 s 69® 
h 45 m 04.3 s -25® 
h 19 m 21.5 s 42® 
h 24 rti 14.5 s 33® 
h 40 m 05.3 s -00® 
h 15 m 00.7 s 38® 
h 39 ra 43.7 s 32® 
h 14 m 36.8 s 56® 
h 33 m 45.8 s 59® 
h 41 in 54.6 s 67® 


(1950) 
55' 17" 

55' 47" 

54' 47" 

33' 46" 

06' 55" 

21* 13" 

13' 44" 

05* 06" 

48' 50" 

30' 17" 

58* 36" 

56' 33" 


Observations 
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Table 5-2 (cont'd) 


Column 1: 
2 ; 
2 : 


4: 


5; 


Object name 

Dates of observations 

Telescope used in observation; either 5 m Hale 
telescope or 4 m Anglo-Australian telescope 
Standards used for calibration of 1 ram flxox from 
galaxy 

—26 2 

Measured 1 mm flux in 10 W/m Hz; errors 


6 : 


are total uncertainty in flux 
Statistical significance of observations of 
galaxy; see text for more details. 


47 


Figure 5-la: Spectrum of NGC 253 from 1 GHz to 10^ GHz, 

including the 1 nan measurement from this work. 
Fitted curve is explained in text. Radio 
points are from Haynes et al. (1975) . Far 
infrared points are from Telesco and Harper 
(1980) , Rieke and Low (1972) , and Hildebrand 
et al. (1977) . The 1 mm point is indicated by 
a triangle. 
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5 

Figure 5-lb; Spectrum of NGC 1068 from 1 GHz to 10 GHz, 
including the 1 mm 3a upper limit from this 
work. The fitted curve has an emissivity 
index of 2, the dotted line corresponds to 
an emissivity index of 1.5. The radio points 
are from Haynes et al. (1975) and Dent and 
Balonek (1980) . The far infrared points are 
from Hildebrand et al. (1977) and Telesco and 
Harper (1980) . 
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Figure 5-lc: Spectrxun of NGC 4631 from 1 GHz to 10^ GHz, 

including the 1 mm measurement from this work. 
The radio points are from Haynes et al. (1975). 
The 1 mm point is indicated by a triangle. 
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Figure 5-ld: Spectrum of NGC 6946 from 1 GHz to 10^ GHz, 

including the 1 mm measurement from this work. 
The fitted curve is explained in the text. 

The radio points are from Haynes et al. (1975) . 
The far infrared points are from Telesco and 
Harper (1980) . The 1 mm point is indicated by 
a triangle. 
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Figure 5-le: Spectrxim of M 82 from 1 GHz to 10^ GHz, includ- 

ing the 1 mm measxirement from this work. The 
fitted curve is explained in the text. The 
radio points are from Haynes et al. (1975) and 
Dent and Balonek (1980) . The far infrared 
points are from Telesco and Harper (1980) . The 
1 mm point is indicated by a triangle. 
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Figure 5-l£: Spectrum of IC 342 from 1 GHz to iO^ GHz, 

including the 1 mm measurement from this work. 
The fitted curve is explained in the text. The 
radio points are from Haynes et al. (1975) . 

The far infrared points are from Becklin et al. 
(1980) . 
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e(v) = Kv^\ (5.2) 

P will be defined as the eraissivity index for the rest of 
this work. Curves of this form have been fit to the obser- 
vations of M 82, IC 342, NGC 253, NGC 6946, and NGC 1068. 
These curves have been plotted in Figures 5-la, b, c, d, e, 
and f. The relevant curve parameters are given in Table 5-3. 
For the other galaxies in our study no far infrared results 
have been published. 


Table 5-3 

Isothermal Dust Model Parameters 


Object 

(“K) 

P 

M 82 

46 

1.5 

NGC 253 

41 

1.5 

NGC 6946 

43 

1.5 

NGC 1068 

33 

2.0 

IC 342 

32 

1.5 (assumed) 


The observations at wavelengths around 100 um are most 
useful in determining the blackbody temperatures, while the 
longer wavelength observations out to 1 mm are better for 
determining the value of p. Observations at wavelengths 
longer than 1 mm are not very useful for studies of the 
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dust, since the flux becomes so small that detection is very 
difficult. In addition, at these long wavelengths disen- 
tangling the thermal flux from the non-thermal synchrotron 
flux is a thorny problem. 

For the galaxies M 82, NGC 253, and NGC 6946 the best 
fit curves have blackbody temperatures of around 40'*K and 
e(v)av^*~. There is not enough data to accurately determine 
the value of the emissivity index for IC 342. A value of 1.5 
was therefore chosen so that a dust temperature estimate 
could be made. NGC 1068 was fit best by a curve with 

•7 

e(v)av“, implying a fundamental difference between this gal- 
axy and the other galaxies in the study. 

The dust responsible for the 1 mm emission from spiral 
galaxies could exist in three different environments. The 
hottest environment (defined here as Region I) would be in 
molecular clouds with embedded stellar or protostellar 
sources. The dust temperature in these clouds has been cal- 
culated to range between 15®K and 100°K, depending on the 
cloud parameters and the location of the dust in the cloud 
(Leung, 1976). Temperatures of this order have been observed 
in such clouds in our galaxy (Thronson and Harper, 1979) . 

Dust in interstellar space (Region II) would have an inter- 
mediate temperature. The equilibrium temperature of un- 
shielded grains in the interstellar medium of a galaxy de- 
pends scrongly on the grain composition. The equilibrium 
temperatures for our Galaxy have been estimated to range 
from 10°K (silicates, Aannestad, 1975) to 40°K (graphite. 
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Werner and Salpeter, 1969). The coldest environment is in- 
side dense dark molecular clouds (Region III). The grain 
temperatures in such clouds have been estimated to be less 
than or equal to 7®K (Leung, 1975). The thermal emission 
from the dust in these clouds would be very difficult to de- 
tect, since the dust temperature is not much greater than 
the cosmic background. 

By combining the far infrared, submillimeter, and milli- 
meter observations from a galaxy, we can determine the rela- 
tive amounts of dust in Regions I and II. Subtracting the 
spectrvun of an average hot molecular cloud from the spectrum 
of a galaxy will leave the contribution due to the dust in 
the colder environments. The average spectral index between 
340 urn and 1 mm for the molecular clouds measured by West- 
brook et al. (1976) is 3.18. This can be compared with 
the spectral index of our galaxies between 340 um and 1 mm. 
Unfortunately, only one galaxy, NGC 253, has been reported 
to have been detected at 340 \im (Harper and Low, 1973) . 

As a result, the 340 um flux can only be interpolated from 
our curves in Figure 5-1 for the other gala.xies we have 
detected. When this is done, the average value of the spec- 
tral index from 340 uin to 1 mm for M 82, NGC 253 , and NGC 
6946 is 3.29. The close similarity of this value with the 
average value for the molecular clouds implies that almost 
all the 1 mm radiation from spiral galaxies is emitted in 
hot molecular clouds. This means that substantial quantities 
of dust outside of these clouds and at interstellar tempera- 
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tures (around 15°K) are effectively ruled out. Since a dust 
cloud at would emit very little radiation, our results 
do not set such stringent constraints in the quantity of 
dust in cold molecular clouds. 

ii) Grain composition ; If the dust grains are all at 
the same temperature, the observed value of the emissivity in- 
dex puts constraints on the grain composition. Aannestad 
(1975) has measured the absorptive properties of various sili- 
cate grain materials, both with and without ice mantles. Of 
the various materials, only fused quartz, with an emissivity 
index of 1.6, is close to our derived value of 1.5. Other 
silicates or silicates with icy mantles have emissivity in- 
dices of 2 or higher. Phillip (1977) has made a similar 
study of graphite grain emissivities. His results for pure 
graphite grains give an emissivity index of 3 for wavelengths 
longer than 100 pm. 

It should be noted that this problem also exists in 
explaining the results from hot molecular clouds, since they 
also have emissivity indices of around 1.5 (Westbrook et al. , 
1976; Werner et al., 1975; Ade et al., 1974). For these 
clouds the low value of the emissivity index is thought to 
be due to the temperature gradient that exists through the 
clouds. The flux contributions of the colder dust located 
at the outside edges of the clouds brings the emissivity 
index down to a lower value than would be indicated by the 
grain emissivity alone. In order to estimate the amount 
of dust and the temperature variations tiiat would be required, 
let us assume that the dust in hot molecular clouds has an 
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emissivity index of 2.5. This is slightly above the average 
for silicates or core-mantle ice silicates, but is well 
within the range of indices for the different silicate 
materials proposed for interstellar dust. Since the eav^*^ 
curves fit the galaxies so well, let us also assume that 
this expression is exactly true. This means that the spec- 
tral shape of the dust emission is given by: 


. 4.5 


rhv/l^T, 


c - 1 


( 5 . 3 ) 


where T_ is the color temperature of the dust, h is Planck's 
constant, and k is Boltzmann's constant. Now the low fre- 
quency tail of this expression can be approximated by 

rT^ 


F (v) a 


,5.5 


,hv/kT _ ^ 


C(T)dT , 


( 5 . 4 ) 


corresponding to our assumed dust emissivity index of 2.5. 
The function C(T) gives the distribution of dust mass as a 
function of temperature. We can find an upper limit for 
in the following way. If equation (5.3) is plotted on a 
log F(v) - log V plot, the slope at a frequency v will be 
given by: 

4.5 (5*5) 

1 -e-^ 


where x = hv/kT^. Similarly the slope of 
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V 

ghv/kT _ ^ 

on the same plot is given by: 

5.5 ^ . (5.7) 

1 - e~*^ 

If the average hot molecular cloud has a color temperature 
of 40'’K, then equation (5.5) gives a slope of 3.3 at 300 GHz 
(1 mm) , the lowest frequency measurement of the dust. In 
order for equation (5.7) to equal this value for the slope 
at 300 GHz, the temperature must run from 40 down to 8®K 
or less. 

Of course, the actual submillimeter and millimeter 
i observations of hot molecular clouds are not accurate enough 

to set such precise temperature limits. Even assuming a 
\ lower limit to the temperature of 15°K, we can obtain excel- 

lent fits to our results. In order to realisticallv esti- 

I 

2 5 

mate the amount of colder eav dust needed to fit the 
observations of our galaxies, we will assume that the dust 
is onlv at two temperatures, and 15“K. A temperature of 

1 

15 is not as arbitrary as it may seem, since this is 
roughly the lower end of the range of temperatures expectec 
in a hot molecular cloud (Leung, 1976) . We will therefore 
fit the galactic spectra to 

F(v) = B(v, 15°K) + C 2 V“*^ B(v, ) (5.8) 
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where Ci and C 2 are different constants for each galaxy. When 
this is done, the parameters given in Table 5-4 are derived. 
These results indicate that our model requires roughly 20 times 
as much 15 ®K dust as color temperature dust in order to fit the 
spectra. It is important to remember that all of this dust 
would be associated with hot molecular clouds in these galaxies. 


Table 5-4 

Two Temperature Dust Model 


Galaxy 

T^ (»K) 

C 1 /C 2 

NGC 253 

34 

12 

NGC 6946 

35 

25 

M 32 

82 

24 


For simplicity, in our derivations of dust masses we will 
assume all of the dust is at an average temperature of 20®K in 
the spiral galaxies. If the actual dust temperatures are 
closer to the color temperatures given in the isothermal mod- 
el earlier, it will only lower our dust mass estimate by a 
factor of 2. Since the value of p is higher for NGC 1068, we 
will use the color temperature for this galaxy. The reader 
should keep in mind that the dust masses that will be derived 
will actually be lower limits, since very low temperature 
dust is being underrepresented. 

Dust Masses 

One mm observations of galaxies are particularly useful 
in estimating dust masses. We have seen that the 1 mm 
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emission is due to thermal radiation from dust at 15-40 “K. 
Since 1 mm wavelengths are on the Rayleigh-Jeans tail of the 
blackbody curve # the expression for the flux received from 
an extended source becomes: 

F(v) = B(v, T(v) (5.9) 

= 2k (vVc^)n t(v) (5.10) 

where B(v, Tjj) is the blackbody function, the dust tem- 
perature, the effective beam size, and t(v) the optical 
depth. 

For a collection of spherical dust grains of number den- 
sity n the optical depth can be expressed as 

t(v) = /mra^Q(v)d£. (5.11) 

In this expression a is the grain radius, Q(v) is the absorp- 
tion efficiency, and d-c is an increment of path length. For 
grains of density p 


T(v) = D ^ (5.12) 

where D is the dust mass column density. 

For an optically thin, isothermal dust region the flux 
from any part of the region can be expressed as 

F(9,i^) = An(9,4>)dw (5.13) 

where A is a constant, and n(8,(^i) is the number column den- 
sity at that position (see Figure 5-2). The signal received 




Figure 5-2: Diagram used for derivation of 1 mm flux. 
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by the telescope is given by 

Fr a /A'n(0,(|))R(0,(J))dfl (5.14) 

where A* is a constant and R(6,4i) is the telescope team 
response. We can define an average column density from 


Fr 

A' /R(d,(p)dSl 


(5.15) 


Now for a Gaussian beam with a FWHM of 55", 

= /R(0,(|))dn a 8 X 10'^ ster, (5.16) 

and corresponds to an effective beam diameter of 66". Com- 
bining all of these equations we find 


_ _ 2 Fr (V) ap 

° “ 3k Q(v) ' (5.17) 

defining D as the beam averaged dust mass column density. 

The value of the factor Q/ap is strongly dependent on 
the grain composition. The calculated and measured values 
for X a 1 mm range from 6 cm /g ("astronom.ical silicate," 
Draine, 1980) to 10”^ cm^/g (graphite, Phillip, 1977). A 
list of some materials that have been suggested for grain 
composition together with their values of Q/ap and p is 
given in Table 5-5. 

This uncertainty in the value of Q/ap is the largest 
source of uncertainty in dust mass estimates from millimeter 
flux measurements. Until the composition of the dust is 
better known, this problem will always be present. For the 


rest of this work we will assume an average silicate value 
of 0.05 for Q/ap at 1 mm. 
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Table 5-5 


Values of Q/ap at 1 mm for Some Suggested Grain Materials 


Material 


(10“^ cmVg) 


Grain Emissivity 
Index / p 


Olivine 

Fused quartz 

Fused quartz with 
ice mauitles 

Lunar rocks 

Pure graphite 

"Astronomical Silicate" 


9.3 2.2 

190 1.6 

48 3 

76 1.9 

10 3 

6000 2 


(From Aemnes tad# 1975; Phillip, 1977; and Draine, 1980.) 


With this value the dust mass column density can be cal- 
culated for the galaxies using their measured fluxes together 
with the temperatures derived earlier. For our 66" effec- 
tive beam size on the 5 m Hale telescope equation (5.17) 
becomes 

D(g/cm^) » (5.18) 

Although we will rather loosely refer to D as the mass col- 
umn density, it is important to remember that this quantity 
has been averaged over our beam and should not be compared 
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directly with other studies using beams of markedly different 
sizes* The results of the dust mass calculations are given 
in Table 5-6. 

Translating these dust mass column densities into total 
masses requires knowing the distances to the galaxies. The 
distances used for these calculations are from various 
sources and are given in Table 5-7. The total mass of dust 
respo'.sible for the 1 mm flux in our 66" effective beam is 
then 

= D(g/cm2)^2(M^^) 3,83 ^ lo® (5.19) 

In this equation D is the dust column density and i is the 
distance to the galaxy. The mass totals are given in Table 

5-8. 

Flux Correlations 

The four galaxies, M 82, NGC 253, NGC 4631, and NGC 6946, 
that were detected in our observations do not form a large 
enough sample to draw statistically significant conclusions 
about the 1 mm flux from spiral galaxies as a class. It 
should also be noted that there was a selection effect in 
choosing which galaxies to observe. The prime criterion 
was the size of the galaxy. In practice, this meant that 
the galaxies were all at roughly the same distance. The 
selection effect entered when the observed galaxies were 
chosen out of this group, since they were chosen because 
their brightness at other wavelengths suggested that they 
might be bright at 1 ram. Thus some of the galaxies may be 
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Table 5-6 

Dust Mass Column Densities 


Object 

Dust Vamperature 

M 82 

20 

1' NE 

20 

1' SW 

20 

NGC 253 

20 

NGC 891 

20 

NGC 925 

20 

NGC 1068 

[33] 

NGC 4244 

20 

NGC 4631 

20 

NGC 5907 

20 

NGC 6946 

20 

IC 342 

20 


(“K) D (X 10“^ g/cm^) 

11.5 

3c < 7.6 
3c < 9.7 
12.9 

3c < 4.4 
3o < 4.6 
3c < 3.5 
3c < 5.8 

4.2 

3o < 7.4 

1.3 

3c < 5.9 


Due to the numerous assxamptions that went into calculating 
the value of D, no formal errors are included. We have 
assumed an average dust temperature of 20 ®K for all of the 
galaxies except NGC 1068 (see text for details) . 

[ ]: color temperature used for this object. 
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Table 5-7 

Galactic Distances 


Object 

Hubble type 

Distance 

(Mpc) 

Ref. 

Diameter 
66” beam ( 

M 82 

Ir II 

3.0 

1 

1.0 

NGC 253 

Sc 

3.4 

2 

(1.4) 

NGC 891 

Sb 

O 

• 

O 

5 

3.3 

NGC 925 

Sc/SBc 

14.8 

6 

4.8 

NGC 1068 

Seyfert 

18.1 

2 

5.8 

NGC 4244 

Scd 

3.8 

6 

1.2 

NGC 4631 

Sc 

5.2 

4 

1.7 

NGC 5907 

Sc 

4.6 

6 

1.5 

NGC 6946 

Scd 

4.2 

6 

1.3 

IC 342 

Scd 

4.5 

3 

1.4 


References ; 

1) Taininann and Sandage, 1968 

2) Sandage and Taitimann, 1975 

3) Baker et al., 1977 

4) Kruinm and Salpeter, 1979 

5) Sancisi and Allen, 1979 

6) Roberts, 1969 

{ ): Diameter of 88” effective beam of Anglo-Australian 


telescope 
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Table 5-8 
Total Dust Masses 


Object 
M 82 
1' NE 
1' SW 
NGC 253 
NGC 891 
NGC 925 
NGC 1068 
NGC 4244 
NGC 4631 
NGC 5907 
NGC 6946 
IC 342 


Total dust mass in 

66" beam (10^ M ) 
o 

5 

3a < 2 
3a < 5 
8 

3a < 16 
3a < 40 
3a < 40 
3a < 4 
4 

3a < 6 

0.8 

3a < 5 
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unrepresentative of spiral galaxies in general. 

In spite of these disclaimers/ it is impossible to re- 
sist the temptation of trying to find correlations between 
our 1 mm results and the results at other wavelengths. Some 
conclusions we can draw immediately. 

1) The 1 mm radiation in spiral galaxies originates 
in thermal radiation from dust. 

2) The 1 mm points are on the Rayleigh- Jeans tail of 
the dust blackbody curves. 

3) The 1 mm, submilliraeter , and far infrared measure- 
ments yield a similar spectral shape for all the galaxies 
detected at 1 mm. This same spectral shape is also observed 
in measurements of galactic molecular clouds. 

As a result of point number 3, there is close correlation 
between the 100 pm and 1 mm flux densities. Our results 
show that for a 1' beam centered on the nucleus of a galaxy 

the average value of is 470. 

^ 100 pm 1 mm 

We can also attempt to find a similar correlation be- 
tween the 1 mm and gas observations. The correlations 
with CO and HI measurements will be discussed later in the 
sections on the individual galaxies. An additional compari- 
son can be made between the 1 mm emission and the total HI 
emission from a galaxy. Such a correlation would imply that 
all spiral galaxies are basically the same, with the only 
differences due to distance and an overall size scale fac- 
tor. There are two ways to check for this. If the source 
of the 1 mm emissiop in galactic nuclei is much smaller than 
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the beam size then the ratio 


i mm 


M. 


HI 


(5.20) 


should be a constant. In this expression F, is the 1 mm 
flux, Myj the total mass of HI in the galaxy, and t the dis- 
tance to the galaxy. If the source of the 1 mm radiation 
is larger than the beam size, then mm/^l should be a 
constant. Using data on the total mass of HI from Roberts 
(1969) , together with the distance estimates listed in 
Table 5-7, it is possible to make these comparisons. The 
results are shown in Figure 5-3. As can be seen in the fig- 
ure, there seems to be no correlation between these quantities. 

One final comparison can be made between the synchrotron 
emission from the galaxies and their 1 mm flux. This com- 
parison is shown in Figure 5-4 and indicates that there may 
indeed be a correlation between these two quantities. If 
this is real, this would imply that there is some physical 
connection between the amount of dust and the strength of 
the synchrotron radiation. This would be understandable if 
the 1 mm emission originates in hot molecular clouds. These 
clouds are the sites of star formation and their numbers 
indicate the general level of activity in the nucleus. A 
high level of star formation would also lead to more super- 
novae which would contribute to the overall level of syn- 
chrotron radiation from the galaxy. Future 1 mm extragalactic 
observations should help to firmly establish this relation- 
ship. 


Figure 5-3: For the galaxies observed at 1 rain the flux is 

plotted as a function of a) the total flux of 
HI divided by the square of the distance and 
b) the total HI mass of the galaxy only. 







1mm Flux CJy) 
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Figure 5-4: The 1 nun flux plotted against the radio con- 


tiruum flux at 2.695 GHz. 
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Individual Objects 
M 82 

The 1 mm flux from M 82 was measured at three points 
along the galactic plane; one point centered on the nucleus 
and the others 1' along the plane on each side of the nucleus. 
The nuclear center given in Table 5-1 is the position of the 
peaJc radiation at 10 ym. Only the flux from the nuclear 
point was strong enough to be detected. The flux measured 
at this point, 1.98 + 0.34 Jy, is in excellent agreement 
with previous 1 mm measurements (Elias et al., 1979 (2.7 + 

0.7 Jy) ) . 

i) Spectrum ; The nuclear region of M 82 is very 
bright in the infrared due to thermal radiation from dust. 

This dust is thought to be heated by large numbers of 
luminous early type stars (Telesco cind Harper, 1980) . 

Since the far infrared spectrum closely matches the spectra 
of galactic molecular clouds with hot embedded sources, the 
infrared radiation from M 82 probably arises from similar 
clouds. As was noted earlier, the 1 mm emission from M 82 
can also be explained as originating in molecular clouds. 

This implies there is very little Region II dust (T^ - 15 °K) 
in the galaxy. 

ii) Gas measurements ; M 82 has been observed in CO by 
Rickard et al. (1973) and in HI by Gottesman and Weliachew 
(1977). The beam sizes for these measurements were 65" and 


3.5', respectively. 
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The CO beam size of Rickard et al. (1977) was very close 
to our 1 mm beam size/ making comparisons between the two 
measurements very easy. In deriving their total CO column 
density Rickard et al. assumed a gas temperature of 26 ®K. 

We have recomputed their results using our dust derived tem- 
perature of 20°K. With this modification their total CO 

i 

column density in the direction of the nucleus becomes 1.5 x 
10^^/cm^. Using their value of lCO)/(H 2 ] = 4 x lO"^ for ' 

the ratio of CO to molecular hydrogen, the column density of 
H 2 becomes 4 x lO'^'^/cm or 1.3 g/cm . | 

The beam size of the HI measurements of Gottesman and | 

I 

Weliachew (1977) was too big for direct comparison with our , 

1 mm results. At this beam size there is no evidence of the 
HI hole found at the nuclear positions in many other galaxies | 

(Krumm and Salpeter, 1979) , although the peak HI emission is , j 

j 

offset 30" to the south from the 10 um peak. The HI peak \ 

corresponds to 1.5 x 10^^ atoms/cm^ or a mass column density 

of 9 X 10~^ g/cm^. This is a factor of 140 lower than the | 

H 2 mass column density estimated from the CO observations 

above. It is interesting that this ratio is so high, since i 

estimates of the inner 500 pc of our galactic center give , i 

i 

[H2]/2[H] = 30 (Rickard et al., 1977). The cause of this ,, 

discrepancy probably lies in the different beam sizes used 
in the measurements. Thus the total gas mass column density 
is close to the H 2 value. 


Combining our 1 mm results with the gas measurements 
gives an approximate gas:dust mass ratio of 110:1. Although 
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this ratio is roughly equivalent to estimates for our galaxy, 
it must be emphasised that the value of (pa/Q) used in cal- 
culating the dust mass column density was only a crude 
estimate. 

iii) Scan comparisons ; Since the CO beam size is so 
close to the 1 mm beam size, the CO and 1 mm seems can be 
compared directly. This is shown in Figure 5-5, The FWHM 
of the CO peak is 80". Since this is only slightly larger 
than the beam size used in the CO study, the nucleus is not 
really resolved. In our galaxy there is a strong correla- 
tion between 1 mm and CO emission. This correlation is also 
evident in our result from M 82; both the CO and 1 mm results 
show strong peaks in the emission at the location of the 
nucleus . 

It is interesting ‘•.o compare these results with the 
results from our own Galaxy. The Galactic center region 
has been mapped at 540 ym by Hildebrand et al. (1978) and in 
CO by Solomon et al. (1972). These results give 540 ym and 
CO source sizes of 13' FWHM (40 pc at 10 kpc) . A similar 
size source in M 82 would be far too small to be resolved 
with our beam. 

In M 82 the distance corresponding to the effective 
bean diameter (1.0 kpc) is only slightly larger than the size 
of the active central region defined by Soliger et al. (1977). 
The non-thermal radio emission from this region is similar 
to that from the Milky Way's galactic center region, although 
the central region of M 82 has a volume roughly 50 times 
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Figure 5- 


> 

I 




: Plot of scans along the plane of M 82 at vari- 

ous wavelengths. The center of the plot is 
the nuclear position listed ‘in Tadale 5-1. The 
1 mm flux is expressed in Jy. The CO results 
are from Rickard et al. (1977) . The 58 u® far 
infrared results are from Telesco and Harper 
(1980) and have a peadc intensity of 1066 ± 44 
Jy. The exact position of the 58 u® peak is 
not known, and the placement of the infrared 
scan is designed to show the size of the fea- 
ture, rather than its exact location. 
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larger than in our galaxy (Soliger et al., 1977). If we 
accept the argument of Soliger et al. (1977) that the nucleus 
of M 82 is similar to the Galactic nucleus » only larger, we 
can estimate hov; large the 1 mm source should be. Since 
Soliger et al. concluded that the linear dimensions of M 82 
are roughly 4 times those in our galaxy, this would correspond 
to a 1 mm source size of 10” at 3 Mpc. 

We can give a definite value for the minimum 1 mm source 
size from 


or 


TT By 2kT^ ^ 
/2 




(5.21) 


(5.22) 


Telesco and Harper (1980) claim to have resolved a 37” 
source of 58 v® radiation in the nucleus of M 82. The 1 mm 
radiation could perhaps originate from the same source. 

iv) Models for M 82 ; What kind of galaxy is M 82 any- 
way? Based on its appearance on photographic plates, M 82 
is classified as an Irll galaxy. The earliest observations 
in the visible and in showed a strange filamentary and 
knotted nucleus. Originally this appearance was explained 
as being due to a recent explosion in the nucleus. More 
recently, observations have uncovered evidence that dispute 
this theory. As a result, other models have been proposed 
to explain M 82 's unusual appearance. Many of these models 
have been summarized by Soliger et al. (1977) and two are 


described below. 
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Gottesman and Weliachew's (1977) work in HI has shown a 
dynamical correlation between M 81 and M 82. They therefore 
proposed that a recent tidal encounter with M 81 has driven 
gas deeply into the nucleus of M 82. This gas caused density 
enhancements that resulted in a massive amount of star forma- 
tion which results in large amounts of infrared emission. 

A somewhat different model has been proposed by Elvius 
(1972) and more recently by Soliger et al. (1977). In these 
models M 82 is moving through a thin cloud of dust and gas 
associated with the M 81 system. This material is sv/ept up 
and accreted into the core of M 82, causing large amounts of 
star formation and dust emission in the nucleus. 

In view of these models it is perhaps surprising that 
the 1 mm results are so normal, with the estimated gas to 
dust ratio and far infrared spectrum similar to our Galaxy 
and other normal spirals. Our work shows that the large 
luminosity of M 82 is due to larger numbers of the same 
type of sources observed in other galaxies. As a result, 
it seems unlikely that there has been a recent massive 
explosion in the nucleus of M 82. Our instrument is not 
sensitive enough to distinguish between the two models 
above. Had it been more sensitive, the thin dust halo pos- 
tulated in the models of Soliger et al. (1977) and Elvius 
(1572) could possibly be detected at 1 mm. 

NGC 253 

NGC 253 is classified as a Sc spiral (Sandage, 1961) , 
showing prominent dust lanes in visible photographs. It is 
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also very bright in the infrared. The 1 mm flux measurements 
from NGC 253 are the only measur«ments presented in this work 
that were made on the 4 m Anglo-Australian telescope. As a 
result of this telescope's smaller size* the beam size for 
these measvurements was 73" FHHM with an effective beam size 
on a uniform extended source of 88". The position of the 
beam center coordinates for the nucleus of NGC 253 were 
approximately 10" away frcxa the 2.2 ym peak reported by Beck- 
lin et al. (1973). 

i) Spectrum: The spectrum of NGC 253 including the 1 

mm point is shown in Fig^lre 5-la. Although the single tem- 
perature fit is not as good for this object as it was for 

M 82* it is obvioi ^ that there is not a great deal of dust 
in the 10"K-30"K temperature range outside of hot molecular 
clouds. 

ii) Gas measurements ; Measurements of the CO in NGC 
253 were made by Rickard et al. (1977) amd Solomon and De 
Zafra (1975) . If we assume a gas temperature of 20^K* the 
results of Rickard et al. (1977) yield a total CO column 
density of 1.3 x 10^^ cm"^. Using (CO]/(H 2 l = 4 x 10”^ we 
find a value of 3.2 x 10^^/cm^ or 1.1 g/cm^ for the nuclear 
column density of H 2 . 

Aperture synthesis measurements of NGC 253 in HI have 
been made by Combes et al. (1977). These measurements were 
made with a 2 ' x 3 ' beam. Their results show a peak HI col- 
umn density of 1.3 x 10^^/cm^ in the direction of the 
nucleus* or roughly 1/3 0th the amount of H 2 derived above. 
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As was the case with M 82, there was no sign of a HI hole 
at the nuclear position. 

Combining the H 2 estimates with the HI and 1 mm observa- 
tions gives a rough gas to dust ratio of 85:1 for NGC 253. 

A short observing schedule precluded euiy 1 mm mapping of NGC 
253. The mapping done in CO by Rickard et al. (1977) shows a 
more extended central region than was observed in M 82. An 
interesting project for the future would be to see if this 
extension exists at 1 mm as well. 

NGC 6946 

Although NGC 6946 is almost face on, it was chosen for 
this study on the basis of the strong central concentration 
observed at other wavelengths (Morris and Lo, 1978; Rogstad 
et al., 1973). The galaxy is classified as a Scd (Morris 
and Lo, 1978). The beam center position for the 1 mm obser- 
vations was approximately 17" off from the nuclear position 
listed in Condon and Oressel (1976). This discrepancy is 
due to an error in the SAG catalogue. The position used in 
our measurements is accurate. 

i) Spectrum ; Although the far infrared data on NGC 6946 

is scanty, there is enough avail^d 3 le to fit the curve shown 

1 5 

in Figure 5-ld. Assuming the emissivity varies as v * , 
this yields a dust color temperature of 43"K. 

ii) Gas measurements ; NGC 6946 has been observed in 

CO emission by Morris and Lo (1978) . Their observations show 
a peak total CO column density of 1.4 x 10^®/cm^ in the 
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direction of the nucleus for an assumed gas temperature of 
20*K. For the purposes of comparison we have reduced this 
data using the technique described in Rickard et al. (1977) . 
For [C0]/[H2] ■ 4 X 10"^ the column density of H 2 is 3.7 x 
10^^/cm^ or 0.12 g/cm^. 

The HI distribution in NGC 6946 has been studied by 
Rogstad et al. (1973) with a beam size of 2*. Their work 
shows a definite HI hole at the position of the visible nu> 
cleus. Undoubtedly this is due to most of the hydrogen 
being in molecular form in the nuclear regions, as implied 
by the sharp nuclear peak seen in CO. As a result, the HI 
column density is only 1.4 x 10 /cm in the direction of the 
nucleus. 

Combining the hydrogen estimates with our 1 mm data we 
derive a rough gas to dust ratio of 92:1 for NGC 6946. Since 
NGC 6946 is almost face on, there is no contribution to the 
1 mm flux from the spiral arms. If the spiral arms have a 
radically different gas:dust ratio than the nucleus, this 
would show up as a difference in the gas; dust ratio for this 
galaxy compared to the edge on galaxies. Such a difference 
is not indicated by our results. It should be noted that the 
100:1 gas to dust ratio quoted for our Galaxy is derived 
from local measurements in the spiral arms. 

NGC 4631 

NGC 4631 is a nearly edge on Sc spiral galaxy. Of the 
14 edge on spiral galaxies studied by Krumm £uid Salpeter 


89 


(1979), N(jC 4631 had tha graatest column dansity of HI. Tha 
distanca to this galaxy runs from 5.2 Mpc to 9.8 Mpc, dapand- 
Ing on the distance determination method used (Krumm and 
Salpeter, 1979) . For the analysis of our worlc we have used 
the lower figure. 

i) Spectrtan : The far infrared radiation from NGC 4631 

has not been detected as yet. It is therefore impossible to 
determine the temperature or emissivity index of the dust. 

We therefore assumed that the dust in NGC 4631 had a tempera* 
ture of 20 *K. 

ii) measurements : NGC 4631 was observed in CO by 

Rickard et al. (1977) but no significant emission was 
detected. Their RMS noise for this object corresponds to a 
3a upper limit on the total CO column density of 1.6 x 10^^/ 
cm**, assuming a gas temperature of 20®K. For (COJ/(H 2 ] * 

4 X 10“^ this meeuis that the 3a upper limit on the molecular 

hydrogen column density is 4.0 x lO^Vcm^ or 0.13 g/cm^. 

The column density of HI has been napped in NGC 4361 by 

Krumm and Salpeter (1979) along the plane of the galaxy with 

a beam size of 3.2*. They find no sign of a hole at the 

position of the nucleus and derive a HI column density of 
21 2 

3.2 X 10 /cm at this position. This gives a lower limit 

—3 2 

to the total gas mass column density of 5.4 x 10 g/cm . 

The r^mge in the gas to dust ratio, going from the 
lower gas limit derived from the HI measurements to the 3a 
upper limit from the CO results, runs from 3:1 to 31:1. This 
indicates that NGC 4631 is an extremely dusty galaxy, or is 
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unusual In some other way that would increase the 1 mm flux 
to CO flux ratio. 

IC 342 

IC 342 is an almost face on Sod spiral galaxy. It was 
chosen to be included in this study due to the data that al- 
ready exists at other wavelengths; data that indicated it 
should be bright at 1 mm. 

i) Spectrum ; Although IC 342 is relatively faint at 

far infrared wavelengths, some broadband detections were ob- 
tained by Becklin et al. (1980). These results, shown in 
Figure 5-lf, are insufficient to determine both the tempera- \ 

ture and emissivity index of the dust. Accordingly, the 

1 5 

emissivity was assumed to vary as v in a meuiner similar 
to the other observed spiral galaxies. Using this assumption 
the temperature is then calculated to be 32**K for isothermal 
dust. The 3a upper limit at 1 mm fits this model. 

ii) Gas measurements ; IC 342 has been observed in CO 
by Morris and Lo (1978) with a 65" beam. Their results show 
a peak total CO column density of 1.1 x 10^^/cm^ assuming a 
gas temperature of 20 "K. These workers have also made a 
scan across the galaxy through the nucleus. As expected, 
the CO emission is strongly peaked at the position of the 
nucleus. Assuming [C0]/[H2l = 4 x 10"^, the column density 
of H2 is 2.8 X lO^Vcm^ or 0.10 g/cm^. 

IC 342 has also been observed in HI by Rogstad et al. 

(1973) using aperture synthesis techniques. Their beam size 
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for the observations was 4'. As was the case with NGC 6946, 
the results show a pronounced HI hole at the location of the 
nucleus. Again, this is probably due to most of the hydrogen 
being in molecular form in this region. The HI column den- 
sity at the nucleus is 7.2 x 10^®/cm^, or l/50th of the 
amount of H2> As a result, our 3a lower limit for the overall 
gas:dust ratio is 17:1 for the nucleus of IC 342. 

NGC 1068 

Seyfert galaxies are characterized by extremely bright 
nuclei. It has already been mentioned that these galaxies 
have been classified into two classes, with a few galaxies 
belonging to an intermediate class (Khachiki^m and Weedman, 
1974; Osterbrock and Koski, 1976). Although this classi- 
fication is made on the basis of observed line widths in the 
visible, it turns out that there is also a difference in the 
infrared spectra between the two classes. NGC 1068 is classi- 
fied as an archetypical Class II Seyfert. Various workers 
have concluded that the infrared radiation from Class I Sey- 
ferts is due to the seune non-thermal source that is responsi- 
ble for the short wavelength radio emission, while the 
infrared radiation from Class II Seyferts is due to thermal 
emission from dust (Neugebauer et al., 1976; Stein and Weed- 
man, 1976). Neugebauer et al. (1976) have carried their con- 
clusions further, postulating that there is a non-thermal 
source in all Seyfert gal2ucies. In their model the non- 
thermal source is stronger in the Class I galaxies than it is 
in the Class II' s; strong enough to vaporize dust grains. 
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Therefore/ infrared emisaion from dust will only be observed 
in Class II galaxies. 

i) Spectrum ! The far infrared and submillimeter spec- 
trum of NGC 1068 cannot be fit with the e(v}av^*^/ isothermal 
dust model that fits the other galaxies observed (Figure S-lb). 
Telesco and Harper (1980) have found a peculiar double hump 

in the far infrared spectrum that can be fit by a two temper- 
ature dust model. The galeucy has not been detected at 
wavelengths longer than 380 ym in the submillimeter. Our 1 mm 
measurement gives a 3a upper limit of 0.99 Jy, over a factor 
of two lower than previous 1 mm upper limits (Elias et al./ 
1978) . 

If the 1 mm emission is due to the lowest temperature 
dust reported in Telesco and Harper (1980) (dust temperature 
» 33”K) / and the dust emissivity is described by e(v)av^z 
then p ^ 2. The very low 1 mm flux also means that there can 
be very little dust in a Region II environment. 

ii) Gas measurements : NGC 1068 has been observed in 

CO by Rickard et al. (1977) . Their CO total column density 

18 2 

for a 33*K gas temperature is 8.4 x 10 /cm . Using 

ICOl/tHj] ■ 4 X 10"^ (which may be greatly in error tor 

22 2 

Seyferts) the column density of H 2 is 2.1 x 10 /cm or 0.70 

2 

g/cm . This column density poses a problem (cf. Rickard et 
al., 1977) in that it implies too much mass within the beam. 
The total mass of stars and gas estimated by dynamical means 
is smaller than the H 2 mass derived above (Walker, 1968). 

The most probable reason for this discrepancy is that the 












value of [C01/[H2] la larger in NGC 1068 than It la in 
apiral galaxiea. 

HI meaaurementa of NGC 1068 have been made by Rogatad 
et al. (1967) and Heckman et al* (1978). Since neither of 
theae obaervationa were made uaing aperture ayntheaia tech- 
niquea, fine acale mapa are not available. With the 3.3* 
beam of Heckman et al. (1978) the average HI column denaity 
in the direction of the nucleua ia only 5 x 10*^ g/cm^. Thia 
value ia ao low aa to be comparable to the maaa column den- 
sity of CO measured by Rickard et al. (1977) . Undoubtedly 
this is due to the different betun sizes used in the two 
observations, the HI flux being much more diluted. As in the 
spiral gal«ucies, we will assume that the HI mass in the 
nucleus is much smaller than the H 2 mass. 

Using the CO measurements only, the 3a lower limit to 
the gas:dust ratio is 200:1. This is very high. Changing 
[C 0 ]/[H 2 ) to a value such that there is no total mass problem 
would bring this ratio down closer to the 100:1 ratio esti- 
mated for our galaxy. 

iii) Models for NGC 1068 : Hildebrand et al. (1977) 

have reviewed the various theories that have been advanced 
to explain the far infrared emission from NGC 1068 and have 
concluded that only thermal radiation from dust can explain 
their observations. Our 1 mm results present a difficulty 
with this conclusion. Since the 3o upper limit at 1 mm is 
so low, the emissivity index of the dust must be greater 
than or equal to 2. This is significantly higher them the 
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«mls«lvlty Indttx o£ other normal spiral galaxies. This means 
that either a) the dust is of a different composition in N6C 
1068 than it is in the spiral galaxies, or b) there is a 
different temperature distribution in the dust emission 
regions in the galaxies and in H6C 1068. 

There is no evidence supporting point a) . In fact> there 
is evidence that the dust is of similar composition in M 82, 
NGC 253, the Galactic center, and M6C li;68 since there is an 
observed 10 ym absorption feature in all of these objects 
(Rieke and Low, 1975; Kleinman, Gillett, and Wright, 1976; 
Gillett et al. , 1975; Aitken and Jones, 1973; Houck et al. 
1980). 

Point b) is more likely to be correct. It has already 
been stated that most of the 1 mm emission from normal 
spiral galaxies originates in hot molecular clouds. In these 
clouds the energy is supplied by small hot stellar or proto- 
stellar objects embedded in the dust and gas. These clouds 
have been treated theoretically by Werner and Salpeter (1969) 
and Leung (1976), among others. These workers have concluded 
that there is only a mild variation of temperature through- 
out the cloud. Nonetheless, the temperature gradient that 
does exist will have the effect of decreasing the spectral 
index of the submillimeter emission. This means that the 
l^u:ge spectral index for NGC 1068 could be due to dust that 
is all at the same temperature. This could happen if the 
dust in NGC 1068 were spread throughout the energy source. 

If this is the case, it is possible to set a lower limit on 
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th« size of the non-thezmal source* 

Hildebrand et al. (1977) have used their sufamilllmeter 
observations and equation (5*15) to derive a lower limit to 
the diameter of the submillimeter emission region of S”* 

At a distance of 18 Mpc this corresponds to 440 pc or a 
volume of 4.4 X 10^ pc^* Since the infrared luminosity of 
N6C 1068 is 3 X 10^^ L. (Telesco and Harper, 1980 ) , the 

w 

energy generation per unit volume is less than lO**^^ ergs/ 
cm • 

If the energy source were non-thermal, it is difficult 
to understand how dust could survive within the source* In 
addition, for dust grains with a radius of 1 ym located in a 
region 440 pc in diameter and with a total dust mass equal 
to the 3a upper limit for NGC 1068 given in Table 5-7, the 
lifetime of a relativistic electron would only be about one 
year before it struc)c a dust grain. This makes it unlikely 
that the dust is mingled in with a non-thermal source* 

One further process for creating isothermal dust is 
given in Harwit and Pacini (1975) . They suggest that the flux 
from active galactic nuclei is due to large numbers of HII 
regions and supernovae* There are two ways such a region 
could have isothermal dust. Large numbers of molecular 
clouds, packed sufficiently close together, would have dust 
temperatures that are more isothermal than the dust tempera- 
tures in an isolated molecular cloud. In other vfords, the 
colder dust in the outer regions of a molecular cloud would 
be heated by the flux originating in a nearby cloud. For a 
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typical HIZ region lumiaoaity of 10^ (Hynn-Nilliana and 
Baeklinf 1974) » it would raquira 3 x 10^ such ragiona to 
account for the infrared luainoaity of MGC 1068* An 
average aized molecular cloud ia perhapa 2 paraeca in diane- 
ter* If we aaaume that a packing denaity of 1 BIZ region 
per 10 pc^ ia neceaaary for iaotheraal duat; the total volume 
occupied by the regiona will be 3 x 10^ pc^. Thia volume ia 
a factor of 10 leaa than the minimum volume derived by 
Hildebrand et al. (1977) . Thia objection could be overcome 
if the packing denaity needed for iaothermal duat were leaa 
than the 1 cloud per 10 pc^ eatimated above. 

There ia a second way that the active region propoaed 
by Barwit and Pacini (1975) could have isothermal duat. Zt 
haa been auggeated by Barwit et al. (1972) that the duat 
inaide of the Stromgren aphere around an ioniaing source 
would be illuminated by an essentially constant flux of Lyman 
photons, irrespective of its location. Zf there were enough 
ionizing sources, the nucleus of MGC 1068 could be one large 
Stromgren aphere with isothermal duat. Following the argu- 
ments outlined in Barwit et al. (1972) and using the lower 
limit to the volume set by Hildebrand et al. (1977) , we find 
that this scenario requires that n^ S 3S0/cm^ within the 
active region. Thia value is within the electron density 
upper limits set by Neugebauer et al. (1976) on the basis 
of the observed emission lines. Therefore, this could be 
what is happening in the nucleus of 1068. 
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Thtt Rarwit-Paelni •eanario dots hava ona sarioua eoa^ll- 
cation. Aa thair modal raquiraa larga amount! of molacular 
cloud! and !tar formation# it will alao maan that thara ara 
fraquant auparnovaa in tha nuclaua of HOC 1068. A typical 
aupaxnova ramnant can aputtar away tha ailicata graina in 
300 of gaa bafora it alowa and coola aufficiantly to ba 
harmlaaa (Oraina and Salpatar, 1979) . With graina that ara 
mora fragila than ailicataa# having icy maatlaa for axan^la# 
avan mora duat can ba aputtarad. If tha auparnovaa occur 
fraquantly anough# thair blaat wavaa will daatroy graina 
faatar than thay can ba craatad in molacular clouda. 

Zf wa lat tg ba tha tima it talcaa a cloud to ttam in- 
aida out# ralaaaing ita graina to tha intaratallar madium# 
than thara muat ba laaa than 




X ^ 


(5.23) 


aupamovaa/yaar for tha graina to axiat. In thia aquation 
Mq ia tha total maaa of duat in tha nuclaua of MGC 1068 and 
wa ara aaatiming a gaa to duat ratio of 100 tl.. Wa can aati- 
mata Asm in tha following way. Lat ua aaauma that lOt of 
tha total luminoaity from NGC 1068 comaa from atara that will 
avantually bacoma auparnovaa. Mow in ona cloud turnover 
lifatima L^t^ anargy will ba ganaratad in tha nuclaua of MGC 
1068# whara L, ia tha total obaarvad luminoaity (3 x 10^^ 

L^# Talaaco and Harpar# 1980). A 10 atar# burning 10% 
of ita hydrogen to helium# will generate 10^2 «rga of anargy 
in ita lifatima. Thua tha number of auch atara that ara 


JSj. ' flfaJISf T l 
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n««d«d to account for tha anargy ralaaaad from tha nuclaua 
of NGC 1068 in tina t^ is glvan by: 

(5.24) 

10*^ arga 

If 101 of thata star* bacoma cupamovaa# thia givaa a aupar- 
nova rata of 


10*^ arga 


(S.25) 


or approximataly 3/yaar. If Mq ■ 2 x 10^ (Talaaco and 
Barpar, 1980) than 

t^ < 2 X 10^ yaara (5.26) 

if tha duat craation rataia to axcaad tha daatruction rata. 
NOW in our Galactic diak tha aatinataa of t^ ranga from 10^ 
to 10^ yaara (Salpatar, 1977} Scovilla at al*, 1979). It 
might ba axpactad that tha cloud tumovar rata in tha nucla- 
ua of NGC 1068 would ba on tha lov* aida of thia ranga « 
ainca tha Harwit-Pacini modal raquiraa a much mora coaqpact 
diatribution of gaa and duat than axiata in tha diak of our 
Galaxy. Tharafora« although aupamovaa aputtaring doaa not 
totally rula out tha Harwit-Pacini acanariO/ it cartainly 
muat ba takan into account in any datailad modal. Sinca 
icy mantlaa aputtar away much mora aaaily than do tha aili- 
cata coraa, tha abova analyaia ahowa that auch cora-mantla 
graina cannot axiat in a Harwit-Pacini modal for NGC 1068. 

Finally# thara ia ona othar way that haa baan auggaatad 
to craata iaotharmal duat. It haa baan ahown (Parch# 1977) 
that tha atallar flux in a planatary nabula will craata 
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enough radiation pressure to drive the dust grains away from 
the central star. As the dust moves away from the star it 
drags the gas along with it. Eventually « a hollow shell of 
neeurly isothermal gas and dust forms uound the steur. Possi- 
bly a similar shell exists in NGC 1068, 2 ilthough on a much 
larger scale. 

One difficulty with any sort of isothermal model is 
caused by the CO measurements. The CO flux level from NGC 
1068 is very high. This implies a great deal of molecul^u^ 
hydrogen. Since so little dust is seen, the gas: dust ratio 
for NGC 1068 is much higher than in our g 2 U.axy. The only 
way to bring this ratio down is to have large amounts of 
very cold unobservable dust. 

Opper Limits . 

The spiral galaxies NGC 891, NGC 925, NGC 4244, and NGC 
5907 were not detected at 1 mm. The various parameters of 
these galaxies, and the results of observations at other 
wavelengths are listed in Table 5-9. The 1 mm derived dust 
mass upper limits assume a dust temperature of 20 *K. 

None of these galaxies have been reported to have been 
observed in CO emission so that estimates of H 2 masses are 
difficult to make. Both NGC 891 emd NGC 925 have been 
observed in HI with sufficiently small beam sizes to be 
useful. NGC 891 is an almost edge on Sb spiral galaxy. 
Sancisi and Allen (1979) have observed this galaxy using 
aperture synthesis technic[ues that realise a beam size of 
25" X 37". A scan along the plane of the gala:^ shows a 
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pronounced hole# roughly 1' in size, in the column density 
at the location of the nucleus. A similar hole was observed 
in N6C 925 by Knunm and Salpeter. Their beam size in this 
observation was 3.2'. Comparing these galaxies to IC 342 
and NGC 6946 , these results would imply that there is a 
great deal of moleculau: hydrogen in the nuclei of these 
galaxies. The HI nucleeu: column densities for both of these 
galaxies are shown in Table 5-9. Since the information about 
the gas is so scanty, no attempt was made to compute the 
gas:dust ratio lower limits for these galaxies. 




Table 5-9 





Upper Limits 



Galaxy 

Type 

Dist 2 mce 

(Mpc) 

3o Upper Limit 
for 0 (g/cm2) 
X 10*3 

2 

(atoms/cm^) 

X 10*^ 

NGC 891 

Sb 

10.0 

4.4 

7.9 

NGC 925 

Sc/SBc 

14.8 

4.6 

1.6 

NGC 4244 

Scd 

3.8 

5.8 

— 

NGC 4907 

Sc 

4.6 

7.4 



HI column densities are from Krumm and Salpeter (1979) and 
Sancisi and Allen (1979) . 

Upper limits on the dust mass column density assume a dust 
temperature of 20®K. 


The Milky Way Galaxy 

It is interesting to consider how our galaucy would 
appear if placed at a distance of 4 Mpc. In some ways it is 
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more difficult to observe our galeocy, since we are located 
in the dusty spiral arms. For example, even the Bubble type 
of our galeucy is not known with any degree of confidence, 
although it is thought to be either Sb or Sc (de Vancouleurs 
and Pence, 1978) . 

At a distance of 4 Mpc a 6€" effective beam size would 
correspond to 1.3 kpc. Thus, we should consider only the 
observations of the central 1.3 kpc of the Galaxy. At a 
distance of 10 kpc, this corresponds to maps of 7 1/2® in 
extent. This is a problem, since most studies, in the far 
infrared emd submillimeter in particular, do not extend over 
such a large area of sky. There is also the problem of in~ 
eluding the contributions from the spiral arms; this contri- 
bution would vary as we change our imagined position outside 
of the Gal 20 cy. 

The largest extent far infrared map of the Galaxy was 
obtained by Hoffmeui et al. (1971) at 100 um. Their mapped 
area, 2® x 3.6®, had a 100 ym integrated flux density of 
7.6 X 10^ Jy. If we assumed that: 

1) all the 1 mm radiation would come from this region, 

2) the dust emissivity would vary as v^*^, and 

3) the dust temperature is 40®K, 

then the 1 mm flux would be 1.5 x 10^ Jy at 10 kpc. At a 
distance of 4 Mpc this would be reduced to 0.1 Jy, or lower 
than the flux from any of the detected galaxies. It might 
be expected that the measured flux at our imaginaury posi- 
tion would be somewhat higher, since our beam would take 
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in more of the galaxy than appears in the map of Hoffman 
et al. (1971) . Our galaxy could therefore be similar in 
brightness to HGC 6946 » which had a 1 mm flux of 0.2 Jy at 
a distance of 4.2 

Summeury and Conclusions 

To summarize our results we find: 

1) The 1 mm flux observed from spiral galaxies is due 
to thermal emission from dust. 

2) The far infrared, submillimeter, and millimeter 
measurements of M 82 euid the spiral galeocies give spectra 
that are similar in shape to those from molecular clouds with 
embedded hot sources that have been observed in our Galaxy. 

The galactic spectra can be fit to an isothermal dust model 
with dust temperatures around 40*K and emissivity indices 

of 1.5, or to a two temperature model with a dust emissivity 
index of 2.5. 

3) There is little dust at temperatures 30*X > T^ > lO^K 
outside of hot molecular clouds in any of the detected 
galaxies. 

4) There seems to be no correlation between the total 
HI mass in a galaxy amd the 1 mm flux from that galaxy's 
nucleus. There may be a correlation between the 1 mm flux 
and the synchrotron radiation from spiral galeucies. 

5) In M 82 the 1 mm emission shows a strong peak at 
the nuclear position, similar to the peak observed in CO. 

6) NGC 4631 seems to be an unusually dusty galaxy. 
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7) Tha submillimater spectral shape of HGC 1068 is 
different than the spectral shape of the other galaxies 
observed. This difference could be caused by differences 
in either the composition of the dust or in the temperature 
structure of the submillimeter emission region. The gas to 
dust ratio is also very high in this galaxy. 

8) The average gas: dust ratio for the galaxies was 
roughly the same as the 100:1 ratio estimated for our Galaxy., 
There eu:e, however# three major sources for uncertainty in 
these measurements-. First of all# the estimates of Q/ap 
were admittedly rather arbitrary. Using a lower value for 
this quantity would bring the gas: dust ratio down. In 
addition# the presence of cold dust (T^ < 10**K) would not 

be detected in this study. Finally# there is some dispute 
about the correct value of [COl/[H 2 l. If we use the value 
suggested in Blitz emd Shu (1980) # the gas:dust ratio would 
be a factor of 5 lower than is reported eUaove. 



Appendix A 
The Bolometer 


Theory 

A bolometer Is a radiation detector that works by absorb- 
ing incident radiation and then measuring the temperature 
rise in the device due to the power absorbed. The ideal 
bolometer will absorb radiation of all wavelengths, giving 
it an infinite bandwidth. Although this ideal is not realiz- 
able, bolometers still offer the largest bandwidths of any 
radiation detectors. 

We cem represent a bolometer schematically as shown in 
Figure A-1. Bolometers are composed of three main elements. 
One element absorbs the radiation, one element measures the 
temperature change, and one element thermally connects the 
bolometer with a cold sink. This last element is needed to 
reset the bolometer after it has absorbed radiation so that 
it can be used to meUce another measurement. These three 
elements are called the absorber, thermistor, and thermal 
link, respectively. In some bolometers some of these ele- 
ments do double duty. For example, in some short wavelength 
bolometers the thermistor also acts as the absorber. 

As will be seen later, making many short measurements 
leads to lower noise than making one long measurement does. 
Therefore, bolometers are traditionally run in am AC mode. 
This meams that the incident radiation varies periodically, 
with the bolometer measuring the magnitude of the variations. 
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Optimisation of bolometers is a complicated problem. 
Following the treatment of Nishioka et al. (1978) we define 
the symbols listed in Table A-1. The DC temperature balance 
of the bolometer is given by: 

+ <*•!> 

It can be shown from statistical mechanics that the thermal 
fluctuations in an object vary as the square root of the 
absolute temperature. Therefore, bolometers become less 
noisy auid more sensitive as they are cooled down. The 
additional temperature variations caused by the periodic 
variability in the radiation absorbed by the bolometer are 
given by: 

AT » ^ (1 +w^c2/g2)“^^^ (A. 2) 

Two different types of thermistor elements are now com- 
monly used to measure the ATrms of the bolometer. The resis- 
tance of superconducting metals, carefully maintained on 
the transition edge temperature, is very sensitive to small 
changes in temperature. Thus superconducting transition 
edge thermistors can be very sensitive. However, bolometers 
using this type of thermistor are more complex and less 
reliable than bolometers using semiconducting thermistors 
(Clarke et al., 1977} . 

Semiconductors with appropriate doping levels can also 
be very sensitive temperature measuring devices. The resis- 
tance of these semiconductors can be described by an energy 
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Table A-1 
Symbol Definitions 

Temperature of cold surface (thermal sink) 

Temperature of bolometer 

Thermal conductivity between cold surface and bolometer 
Heat capacity of bolometer 

Total power into bolometer from absorbed radiation 
and resistive heating 

2irX chopping frequency 

RMS temperature change of bolometer 

RMS change in radiative power absorbed in bolometer 

Current through bolometer 

Resistance of bolometer 

Energy gap of thermistor crystal in *K 

Radiation absorbed by bolometer 

Boltzman's constant 

Plemck's constant 

Resistance of load resistor 

Temperature of load resistor 

Excess current noise constant 

Preamp current noise 

Preeunp voltage noise 

Bolometer absorption efficiency 

^ ^B 
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gap relation 

R - Ro«^^^ (A. 3) 

where A depends on the doping level. The bolometer used in 
the Cornell ^e system used a Germanium crystal doped with 
indium and Mtimony as the thermistor element. The resis- 
tance of the thermistor element was measured by connecting 
it to a current source and then measuring the voltage drop 
across the element. The simplest current source is a 
voltage source and a large load resistor connected in paral- 
lel. With this current source the voltage response of the 
bolometer becomes: 


S(volts/watt) 


I dRjj/dTj, (1 - 6) 

Ggd +w2c2%^2j1/2 


(A. 4) 


In this equation 


Ge - G-I^ 5^ (1-26) (A.5) 

This effective thermal conductivity reflects the extra heat- 
ing by the bias current as the resistance changes. Again, 
following the treatment of Nishioka et al. (1978) we can 
write the voltage noise in the output of the bolometer as 
the combination of the various noise sources added in 
quadrature. 


'tot 


s*n|p + 


* 


N 


JB 


+ n: 


EB 


+ Njp + N, 


VP 


(A. 6) 


In this equation, Nsp is ths noise to statistical fluctua- 
tions in the baokground power and 
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'^.^hvAT 


s /iP]i}kT for hv « kT 


(In this equation T is the temperature of the background.) 
NfP is the noise due to statistical thermal fluctuations in 
the bolometer and 


'4kTjj^G 


(A.8) 


Njl is the Johnson noise in the load resistor and 




(A. 9) 


Njg is the Johnson noise in the bolometer and 

Ngg is the excess current noise in the bolometer and 


This noise is also called contact noise and arises from poor 
electrical contacts. Zt is impossible to eliminate, but 
with proper construotion techniques the constant N can be- 
come very small. 

N^p is the excess current noise of the preamp and 


(A.12) 




uo 


and **VP la tha voltaga noise of the preanp and 

■ «n/*^ 

The total power, P, into the bolometer that imast be carried 
away by the thermal link is given by 

P - Pb l^»b 

Looking at these equations a number of general points 
can be made about bolometer optimisation. 

1) It is desirable to have as small a heat capacity 
as possible. In pritotice this means going to low tempera- 
tures, since the heat capacity varies as for Insulators 
and superconductors and as T for conductors. 

2) Even assuming the heat capacity can be made very 
small, the total power, P, into the bolometer deteannines 
the minimum value of G that can be used. From equation A.l 
we see that too small a value of G means that the bolometer 
temperature will rise to unacceptable values. Therefore, 

it is important that the background radiation falling on the 

bolometer be filtered so that the bolometer is heated only 

by radiation wavelengths of astronomical interest. Further- 

2 

more, it is also important to keep the electrical I R 
heating of the bolometer at a minimum. 

3) The maximum value of is limited by the proper- 
ties of the JFETs used in the preamps. For most transis- 
tors of this design the current noise and voltage noise 
will be equal with an input Impedance of approximately 1 M(2. 


Ill 


4) Using ths snsrgy gap rslation of aquation A. 3 for 
tha bolomatar wa find 


dRjj ARjj 


(A.15) 


Th^rafora, kaaping tha bolomatar at as low a tamparatura as 
possibla is is^ortant for high sansitivity. 

5) Tha load rasistor should bs as larga as is physical- 
ly convaniant. 

Zt can ba saan that bolomatar optimisation is a compli- 
catad problam, with paak rasults obtainad only for a givan 
bolomatar anvironmant. Tha bolomatar usad in tha Cornall 
^Ha systam was carafully optimisad for tha background powax 
load« praamp propartias, and cold surface tamparatura of 
tha dewar being used. Tha electrical performance of tha 
bolometer is givan in Table A-3. Tha optical performance of 
tha bolomatar will ba covered in Appendix B. 


Total System Noise 

Some of the various sources of boloaeter noise have 
been discussed above > To gat a conqplata picture of tha sys- 
tam performance it if: useful to consider all of tha noise 
sources* The noise sources can be divided into four main 
categories . 

1) Dewar noise* Dewar noise incorporates all of the 
terms in equation A* 6 as well as additional noise caused by 
microphonics and changing sensitivity* The measured magni- 
tude of the various terms of equation A* 6 is given in Table 
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A-2. In addition, tha dawar optical raaponaivity and 
noiaa ara plottad in Figura A-2 aa a function of bolomatar | 

biaa currant* j 

Microphonics is always a problam with high impaduica 
datactor systams such as bolomatars. In ordar to miniaisa 
this problam, avery wira was sacuraly gluad down insida tha . 

dawar whila tha dawar was constructad to ba as rigid as j 

possibla. Thasa procaduras causad tha rasidual microphonics \ 

to ba of a high rasonant fraquancy (v » 10 Hz) and to ba 
quickly damped* A Fouriar analysis was mada of t);a total 
dawar noisa and tha results ara shown in Figura A-3* It 

i 

can ba seen that there is soma 1/f noisa, but at tha 11 Hz 

chopping frequency its magnitude is small. , j 

■ I i 

In ordar to estimate tha degree tha sensitivity of the ! 

dawar changed during an observing run, careful maasuramants ‘ ! ' 

ware made of tha system response in tha laboratory as a j ^ 

1 1 

function of time and tilt angle of tha dawar* It was found ! 

that tha dawar sensitivity remained constant to within a few |; 

! I 

percent as long as tha dawar was not tipped more than 70* 

r 

from tha vertical and as long as sufficient time was allowed I 

I 

for the dawar to reach operating temperature after tha heat 
switch was closed. In practice, this tharmalisation time 
was only about 20 minutes* | 

2) Photometer noise. Photometer noise arises from two | 

different sources* If the chopper throw is sat longer than 
about 4 1/2', tha wobbling tertiary mirror will begin to 
sti^ke the counterbalanced magneto* This causes a serious 


Table A>2 


Dewar Noise at 10 Hz, Measured at the Bolometer 


Background power noise 
Thermal fluctuation noise 
Bolometer Johnson noise 
Load resistor Johnson noise 
Excess current noise 
Preamp voltage noise 
Preamp current noise 


5.2 nV//Hz 

5.0 nV//Hz 
5.6 nV//Sz 

3.2 nV//Sz 
4.9 nV//gz 

4.1 nV//Hz 
5.4 nV/vfz 


Table A-3 

Measured Bolometer Electrical Performance 


Time constant C/Ge = 13.9 msec 

Heat capacity C * 1.50 x 10”^ J/*K 

Effective thermal conductivity Gg = 1.08 x 10“^ W/*K 

Thermal conductivity G « 0.91 x 10"^ W/®K 

Operating temperature T]^ » 0.43^K 

Load resistor > 36.5 MD 


iP, i W.ertsrr 




< * • 
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Figure A-2: 


I 

I 

i 

Bolometer optical signal and noise as a £unc- • | 

tion o£ bolometer bias current. The optical 

% 

signal was measured using a 10 Hz chopped 
liquid nitrogen source with eui f/35 beam. 



optical 10 Hx 



Signal 

(mV) 


1 






Figure A*>3: Fourier transform spectrvim of noise from 

bolometer and preamp. The high frequency 
components have been filtered so that there 
was no aliasing problem. 


Nois# f Arbitrary Unltal 
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distortion of the mirror position waveform and introduces 
large amounts of microphonic noise in the dewar. This noise 
source ceui be eliminated by setting the chopper throw at 4* 
or less. Another possible source of noise from the photome- 
ter is important only during observations of bright sources. 
If the wobbling mirror is not perfectly periodic in time, 
the beam position on the sky will change. Since the beam is 
so large, this effect is very small, but could contribute to 
the noise from a bright source. For fainter objects the 
other sources of noise are more important. 

3) Atmospheric noise. Due to atmospheric opacity at 
submillimeter and millimeter wavelengths, atmospheric noise 
is significant in all but the very best weather. Sky noise 
is very poorly understood, but is thought to have roughly a 
1/f frequency dependence . During our 1 mm observations the 
sky noise is due to inhomogeneities in the atmospheric water 
vapor. As these inhomogeneities drift through the telescope 
beam they will cause variations in the optical depth of the 
atmosphere at 1 mm. The variability will cause noise from 
both the atmospheric emission and treinsmission. The noise 
caused by the variability in the atmospheric transmission is 
only important for the brighter sources, while the varia- 
bility in the emissivity is more important when observing 
faint sources. 

Sky noise is difficult to parameterise , since its 
effects vary with the observing site, telescope, and weather 
conditions. In general, the lowest sky noise occurs on 
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days with very little water vapor. The worst sky noise in 
a visibly blue sky occurs before a frontal system moves in. 
Any type of water droplet clouds such as cumulus or stratus 
clouds render observations useless. Cirrus clouds, being 
composed of ice crystals, are a problem only in that they 
can obscure the bright stars needed for the offset pointing 
technique used during the observations. In fact, cirrus 
clouds have been present during some of the driest days at 
Mount Palomeu:. With a blue sky, the precipitable water 
vapor levels must be lower than 6 mm for photometric results. 

4) External man-made noise. The other noise that has 
been noted during some observations has been caused by RF 
interference. At Mount Palomar a nearby radar station cem 
cause RF pickup when the dome slit is aimed due east. At the 
Anglo-Austral ian telescope a television broadcast repeater 
located on an adjacent peeUc gave some excess noise. By care- 
fully shielding the system electronics these problems were 
minimised, but were always present in some degree when the 
respective transmitters were on the air. Luckily, neither 
of these transmitters ran continuously, so that low noise 
observations were possible when they shut off. 

Bolometer Construction 

It has been shown that it is desirable for bolometers 
to have the lowest heat capacity possible. The easiest way 
to reduce the heat capacity of the bolometer is to make it 
smaller. However, beyond a certain point further reductions 
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j 
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in size are no longer possible, since the bolometer will 
become less efficient in absorption when its size approaches 
the wavelength of the radiation. 

Traditionally, bolometers for the detection of wave- 
lengths of ap to 100 ym have been constructed of semiconduct- 
ing single crystals with two soldered leads providing support, 
a thermal connection with the cold surface, and electrical 
connections to the semiconductor. The incident radiation 
is absorbed either in the bulk crystal itself or in a thin 
layer of absorbent material painted on the crystal surface. 

For longer wavelengths this type of bolometer is no 
longer suit^U^le for three reasons. 1) The heat capacity 
of semiconductors is fairly high for bolometers large 
enough to absorb efficiently. 2) The bulk crystal becomes 
transparent at these longer wavelengths. 3) Any ed^sorbant 
material painted on the bolometer would require a thicJcness 
that would increase the heat capacity dramatically. There- 
fore, for wavelengths longer than 100 ym composite 
bolometers are more sensitive than bulk semiconductors with 
soldered leads. 

In a ccmiposite bolometer different specialised parts 
are assembled together so that the heat capacity is at a 
minimum while the physical dimensions are large enough for 
efficient absorption. The composite bolometer constructed 
for 1 ram astronomy at Cornell is shown in Figure A-4, while 
Table A-4 lists the sizes and heat capacities of the various 


elements 
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The basic structural element of the bolomter is a 

sapphire crystal chip 4.8 mm x 4.8 mm x 0.13 mm in size. 

Sapphire was chosen because its high Debye theta will cause 

it to have an extremely small heat capacity at the dewar 

operating temperature. Diamond has an even smaller heat 

capacity, but is much more expensive. In any case, even with 

sapphire the other elements in the bolometer provide the 

2 

dominant heat capacity. The 23 mm area of the sapphire is 
large enough to provide efficient absorption of 1 mm wave- 
length radiation. 

A layer of Nichrome has been evaporated on one side of 
the sapphire chip. Nichrome was chosen because its resis- 
tivity is high enough that continuous films of relatively 
high impedemce are easily obtainable, because the resistivity 
changes very little with temperature, and because the alloy 

3 

does not become a superconductor at pumped He temperatures. 
The film thickness was chosen so that the surface resis- 
tance of the metal film matched the impedance of free space 
inside the sapphire. With this criterion there will be no 
multiple internal reflections when the radiation is incident 
on the sapphire surface, travels through the sapphire, and 
is absorbed in the metal film. The absence of these internal 
reflections also means that the adsorption will be indepen- 
dent of frequency. With frequency independent absorption, 

50% of the radiation is absorbed, 25% is reflected, and 25% 
is transmitted on a single pass through the bolometer 
(Clarke et al., 1977). 
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by 


Th« lmp«danc« of fr«« ipaco inalda tha sapphlra la givan 


Z 


n - 1 


(A.16) 


whara Zq la tha impadanca of fraa apaca in a vacuum and haa 
a valua of 377 Q/Q (Clarka at al., 1977). For aapphixa at 
mm wava fraguanclaa n « 3 ao that a aurfaca raalatanca of 
189 Q/Q waa raquirad for tha Michroma film. With tha raaia- 
tivity of Nlchroma of 1.08 x lO”^ Q-cm wa find a film thick* 
naaa of 571 nm ia naadad for tha proper impadanca. In 
practical the thicknaaa of tha Nichroma film waa monitored 
during tha evaporation procaaa by maaauring tha viaibla light 
tranamiaaion through tha film. In tha dawar tha bolometer 
ia mounted in an integrating cavity ao that almoat all tha 
radiation ia eventually eUaaorbad in the bolometer. 

In order to meaaure tha temperature change of tha 
sapphire-Nichroma combination a amall cryatal of germanium 
doped with indium and antimony waa epoxiad onto the aapphire 
surface. Thia cryatal waa supplied by A. Sievers and its 
properties are described in Drew and Sievers (1969) . 

The crystal was prepared in tha following way. First, 
it was cut to approximate size using a wire saw. The rough 
crystal was than lapped on glass plate using a slurry of 
water, glycerine, and #600 grit Unasil abrasive. After being 
washed in acetone and distilled water, the crystal was 
etched in a solution of CF*4 until it was shiny and smooth. 
(CF*4 is composed of a 5:5:3 ratio of 70% nitric acid, 48% 
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hydrofluoric acid, and glacial acatic acid# raapactivaly. 

A f«w drops of 481 hydrobromic acid wars addad to aach 100 ml 
of mixture in order to spaed the etching process.) Whan the 
etching was complatad# the crystal was immediately quenched 
in anhydrous methyl alcohol. 

The first step in the soldering procedure was to wet 
the ends of the germanium crystal with pure (99.99%) indium. 

A soldering iron was plugged into a Variac so that its tem- 
perature could be maintained at just above the melting point 
of indium. The germanium was held with tweezers with a 
very small droplet of Special-X soldering flux on the bottom- 
most end. Touching the germanium briefly to the indium- 
tinned and fluxed soldering iron tip caused the end of the 
crystal to wet with indium. The procedure was repeated with 
the other end of the crystal. At this time the fine wire 
brass leads were also tinned at their very ends with indium. 

The crystal was then placed on a resistance strip 
heater and the fine wires mounted as shown in Figure A-5. 
Small droplets of Special-X flux were also placed in the 
positions indicated. While the process was monitored under 
a microscope# the resistance strip heater was turned on 
just long enough to cause the fine wires to be soldered to 
the ends of the crystal. The vapors from the boiling drop- 
lets of flux provided the necessary fluxing action. 

The completed crystal with leads attached was washed 
thoroughly in anhydrous methyl alcohol to remove all traces 
of the remaining fltix. The germanium crystal was then tested 
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Flgur* A-5: Overhead viaw of 6a thamiator on raaiatanca 

haating atrip juat prior to aoldaring. 
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at ^He temperatures to make sure there was very little excess 
current noise. This noise , which is 1/f in character, arises 
predominantly from poor electrical contacts. After passing 
this test the crystal of germanium was glued into place on 
the sapphire crystal chip using Emmerson and Gumming Stycast 
2850-?? epoxy resin mixed with catalyst #9. 

The fine brass wires, being only 0.0008" in diameter, 
were much too fragile to support the weight of the bolometer. 
Therefore, strands of Dacron 15 ym in diameter (individual 
strands of unw^uced dental floss) were epoxied to the four 
corners of the sapphire chip. The other ends of these 
strands were epoxied to the bolometer mounting ring. Fur- 
ther rigidity was provided by additional Fyrex fibers 
laced around the bolometer euid attached to the mounting 
ring. One of the brass bolometer leads was soldered to the 
mounting ring using Wood's metal to form the ground point 
for the bolometer. The other lead was epoxied to the ring 
with Teflon sheets insulating the wire from the ring. Using 
this mounting system for the bolometer gave enough rigidity 
to reduce the microphonic noise to a tolerable level. The 
mounting ring was fashioned from brass and was cemented to 
the OFHC copper liquid ^He chamber. 


Appendix B 

Filters and Field Optics 

As was noted in Appendix A, it is most important that 
the incoming radiation entering the deweu: be filtered so 
that only the wavelengths of astronomical interest are ab- 
sorbed by the bolometer. As is seen in Figure 2-1/ sub- 
millimeter wavelengths shortward of approximately 750 ym are 
almost totally absorbed by the atmosphere at ground based 
observing sites. As a result/ the filter stack shown in Fig- 
ure B-1 was used to form a low pass filter with a half power 
point of approximately 750 ym. Taking the various elements 
separately: 

1) The Teflon window. The primary purpose of the Teflon 
window is to provide a vacuum tight window that is trans- 
parent to 1 mm radiation. The transmission of Teflon in the 
sxibmil lime ter is given in Hadni (1967) and Muehlner and Weiss 
(1973) . Using these results we find that the transmission 

of the Teflon window is approximately 90% from 750 ym to 
1500 ym. 

2) Black polyethylene. The purpose of the 77*K black 
poly filter is to eliminate the visible and near infrared 
radiation. Hadni (1967) gives the transmission of black poly 
as 80% in the 750 ym to 1500 ym range. It is opaque in the 
visible and near infrared out to 16 ym. 

3) Fluorogold. The 4.2‘’K fluorogold filter provides 
additional far infrared filtering. The transmission curve 
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for 3/32" thick fluorogold rises rather gradually longward 
of 280 um with a half power point of 400 um (Muehlner and 
Weiss ^ 1973). Although this material is opaque from 1 ym to 
40 ym, it is barely translucent in the visible, hence the 
need for the black poly. In the 750 ym to 1500 ym range the 
transmission of the fluorogold filter is 90%. 

4) KRS-5. At room temperature KRS-5 is opaque from 40 
ym to at least 1500 um. However, when it is cooled to 4.2"K 
it opens up at longer wavelengths with a half power point of 
750 ym for a crystal thickness of 1 mm. The biggest diffi- 
culty in using this material as a filter for millimeter work 
is caused by its very large refractive index (n = /e » 5.7, 
Wolfe, 1965) . The large index causes significant reflection 
losses from the filter sxir faces. In an attempt to reduce 
these losses, quarter wave (at 1 mm) plates of Infrasil 
quartz were mounted on each side of the filter. Although In- 
frasil quartz has a refractive index of only 2.0 at 1 mm, it 
was deemed to be the best material available due to its low 
water content aund structural strength. Ideally, the quarter 
wave plates should have had a refractive index of /577 or 
2.39. Even with the antireflection plates, the transmission 
of the KRS-5 filter is only 60% at 1 mm. 

Multiplying the transmission of the veurious components 
of the filter stack we arrive at a total transmission of 
40% at 1 mm. 

The field aperture lies just behind the filter stack. 
This aperture performs two functions. First of all, it 
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determines the beam size on the shy. The 1/4" aperture 
corresponds to 64" on the sky with the 5 m Hale telescope. 
Secondly, the aperture size is roughly the same as the size 
of the Airy disk for an f/4 beam at 1 mm. As a result, the 
aperture acts as a long wavelength cutoff filter. The expres- 
sion for the fraction of light contained inside a radius r 
of the Airy pattern for a circular aperture is given by: 



where k * 2ir/X, a is the radius of the circular aperture 
(100" for the 5 m Hale telescope) , R is equal to where 

D is the distance from the aperture to the Airy pattern (800" 
for the 5 m telescope at f/4) , and Jq and eure Bessel 
functions of order zero and one, respectively. Using r “ 

1/8" and the above values for the Hale telescope, we find: 

P(X) » 1 - Jq2(2.48/X) - J3^2(2.48/X) (B.2) 

where X is measxxred in millimeters. A plot of this function 
is shown in Figure B-2. 

It should be noted that this expression is not exactly 
correct, since it does not take the obscuration due to the 
prime focus cage into accoimt; however, the change resulting 
from introducing this factor is only of the order of a few 
percent. Furthermore, equation B.l is only valid for a 
point source; for extended sources there is no Airy pattern. 
As a result, the long wavelength response of the dewar 


Figure B- 


: Transmission of the various apertures and 

combinations of apertiures in the dewar and 
telescope. 
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changes somewhat depending on the nature of the source. 

A short focal length fused silica field lens is located 
behind the field stop* The focal length was chosen so that 
together with the integrating cavity entrance aperture the 
f/f of the dewar would match the f/4 beam of the telescope- 
photometer combination. The transmission of the fused silica 
at 1 mm is somewhat uncertain, since the water content of 
the material is not known. Assuming that it is similar to 
the fused silica measured in Hadni (1967) a value of 75% can 
be assigned for the lens transmission. 

The long wavelength attenuation due to the integrating 
cavity entrance aperture can be treated in a similar manner 
as the field stop* In this case the values to be used in 
equation B«1 are: a « 3.18 mm, R ■ 5.28 mm, and r « 0.66 mm. 

Since the beam f/#'s are the same, equation B.2 is still 
valid. There are, however, two additional complications: 

1) The aperture is roughly equal in size to the wave- 
length of the incoming radiation. The transmission of nor- 
mally incident radiation through a circular aperture is 
given by: 


21cA 

/o 

where a is the radius of the aperture and k ■ 2:r/X. The 
integrating cavity entrance aperture has a diameter of 1.32 
mm, hence equation B.3 becomes 
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X ,16.59/X 

wh«r« X is in millimeters. This function is plotted in Fig- 
ure B-2 . 

2) Radiation will not illuminate the field aperture 
evenly if a point source is being observed but will instead 
have an illumination pattern described by the Airy func- 

I 

tion. Since the field stop is large enough so that most of | 

the Airy central peak fits inside of it at 1 mm, this is a | 

! 

significant effect. Therefore, the expression given in equa- 
tion B.2 is no longer valid unless an extended source is 
being observed. Unfortvinately, the magnitude of this effect 
cannot be measured in the laboratory with our present test 
apparatus since the equipment illuminates the field aperture 
evenly. Mathematically, the problem is very complex, so 
that it would need to be solved numerically. Observationally, 
many of our sources were neither point like nor greatly ex- 
tended, so that the problem would have to be solved as a 
function of source size as well. Due to these complications, 
we have decided to approximate the illumination of the 
field stop as being totally even, which allowed us to use the 
measured dewar response in the lab. 

Combining the transmission of the filter stack, field 
lens, euid Integrating cavity entrance aperture at 1 mm we 
find a total dewar transmission of 22% at 1 mm. 
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The measured spectral response of the dewar was meas- 
ured on both a grating monochrometer and a Fourier transform 
spectrometer* These results are presented in Figures B-3 
and B-4. It can be seen that the grating monochrometer had 
much better resolution. The mercury vapor source of the 
grating monochrometer was assumed to have a flat spectral 
output in the 500 ym-1500 ym range. In addition, special 
precautions were taken to eliminate the effects of diffrac- 
tion in the monochrometer. The results of the Fourier trans- 
form spectrometez have been divided by the response of a 
Golay cell, which was assumed to have a flat response. 

Now that we have estimated or measured the transmission 
and spectral response of the various optical elements of the 
dewar we are in a position to be able to give the optical 
performance of the bolometer and dewar. The dewar optical 
NEP was measured on a test stand consisting of a variable 
aperture, chopper wheel, and liquid nitrogen source. In 
both the optical and electrical NEP measurements the noise 
was measured at a frequency of 10 Hz. Integrating the dewar 
response curve over the black body function for the liquid 
nitrogen chopping gives the total power entering the bolome- 
ter integrating cavity as 


P 


3.86 X 10“9 

(f7#P 


watts 


(B.5) 


where (f/#) is the f number of the test stai.d aperture 
(dewar-aperture distance/aperture diameter) . This expression 









Figure B-4: Fourier transfora spectroneter dewar response. 
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is valid only up to £/4, the entrance aperture of the dewar. 
Now diffraction has the effect of smoothing the edges of the 
dewar f-beam pattern. Therefore, in order to accurately 
measure the optical performance of the bolometer, f/35 
source was used in the test stand. Finally, in order to com- 
pare an optical NEP with the electrical NEP the optical NEP 
must be converted to a DC response. This can be done using 
equation A. 2. With the time constant given in Table A-3 
and a chopping frequency of 10 Hz this amounts to a factor 
of 0.75 that must be multiplied against the 10 Hz optical NEP 
measurement to convert it to a DC measurement. 

The optical performance of the dewar and bolometer is 
given in Table B-1. This performance cem now be compared 
with the electrical NEP derived in Appendix A to find a 
bolometer eUssorption efficiency of 37%. 


Table B-1 

Bolometer Optical Performance 

Dewar optical NEP (10 Hz) 2.5 x lO”^^ W/vliz 

Bolometer optical NEP (10 Hz) 5.1 x 10"^^ W/vIjF 

Bolometer optical NEP (D.C.) 3.9 x lO"”^^ W//5F 

Bolometer absorption efficiency 


37 % 


Appendix C 

AtiROspheric Corrections 

In order to correctly interpret the 1 mm signals meas- 
ured on the telescope / the effects of the atmosphere must be 
removed. This is done by modeling these effects using an 
appropriate atmospheric model, deriving a number that is a 
measure of the magnitude of the atmospheric effects, atnd 
dividing the signals received on the telescope by this number. 

Water vapor is responsible for most of the atmospheric 
absorption at 1 mm. The water vapor absorption parameters 
used in our atmospheric model were derived from Burch (1968). 
The data correction method that was employed was exactly the 
same as the one described in Elias et al. (1978). 

The flux density from any object is approximated by 

F(v) » Fq(v/vq)“ (C.l) 

where Vq is 3 x 10^^ Hz and n is the spectral index of the 
object. Therefore, the total system response to this object 
will be given by 

S * /P(v)R(v)dv = Fq/ ( v/vQ)*^R(v)dv (C.2) 

where R(v) is the total dewar, telescope, amd atmospheric 
response at frequency v. We can therefore define a function 

Q 5 /(v/vo)"H(v)dv (C.3) 

The function R(v) will depend on the atmospheric 
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response, which is predominantly a function of the water 
vapor in the line of sight, and the dewar and telescope re- 
sponse, which will vary slightly with the spatial extent of 
the source. Therefore, Q(n,w) can be tabulated for both point 
and extended sources, where w is a measure of the water vapor 
in the line of sight. Some representative Q curves are 
shown in Figure C-1. The curves were derived from the dewar 
response shown in Figure B-3. 

During 1 mm observations the water vapor is measured 
hourly with the Westphal hygrometer as long as the svin is 
up. During the night the water vapor is estimated by com- 
paring the observed signals of standeurd objects with earlier 
measurements of the same objects measured in the daytime. 

Since the sensitivity of the bolometer remains constant to a 
high degree, the atmospheric model can then be used backwaxds 
to find the precipitable water vapor. 

The accuracy of the Q curves cam be checked in two ways. 
The results along a curve can be checked by following an 
object as it seta, since this will cheuige the airmass and 
hence the water vapor in the line of sight. A plot of the 
raw signal from NGC 7027 versus water vapor integrated along 
the line of sight is shown in Figtire C-2. Although there is 
considereUsle scatter in the points, there is little evidence 
of a steep systematic slope in the data. There may be a 
very slight negative slope, but the errors introduced by this 
are much less than the other calibration errors. 



Figure C-1: Plot of Q as a function of water vapor, w, for 

the two spectral indices and both extended and 
point sources. 




Figure C-2: Plot of signal from NGC 7027 divided by Q for 

different water vapor levels. 



i 

I 
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The curves can else be checked by observing two objects 
with different spectral Indices and with known fluxes. We 
have used the planets and NGC 7027 for this teat. Interpolat- 
ing from radio and far Infrared measurements we derive a 
1 mm flux from NGC 7027 of 5.5 Jy (Terzlan, 1978; Moseley, 
1979) . Using our method of data reduction and our atmospher- 
ic model we find 2ui observed 1 mm flux of 5.8 ± 0.8 Jy from 
NGC 7027 with Saturn as a calibrator. Since the spectral 
Indices of NGC 7027 and Saturn were taken to be 0 and 2, 
respectively, this agreement between our flux and the Inter- 
polated flux shows that the Q's are also accurate between 
curves. The 1 mm measurement of Elias et al. (1978) gives a 
value of 4.8 ± 1.0 Jy for NGC 7027 at 1 mm, in good agree- 
ment with our results. We therefore conclude that our atmos- 
pheric model does correct for most of the effects of the 
atmosphere . 

Using the Q values derived from our atmospheric model 
we find that the 1 mm flux from am unknown object is given by 



where Fg is the flux from the 1 mm standard, Wg is the water 
vapor along the line of sight to the standard, ng is the 
spectral index of the standard, and Sg is the raw signal 
received from the standard. S, w, auid n are similarly de- 
fined for the object. 
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We can define the effective wavelength during an obser- 
vation from 


^eff - 1 

Xeff plotted for two values of n as a function oi w for 
two values of n in Figure C-3. This figure shows that the 
effective wavelength can be considered to be 1 mm for just 
about any object in any water vapor conditions. 
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3 

Figure C-3: Plot of the effective wavelength of the He 

millimeter system as a function of water 


vapor, w 
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